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Abstract 

We study the antiferromagnetic (/-state Potts model on the square lattice 
for q = 3 and q = 4, using the Wang-Swendsen-Kotecky (WSK) Monte Carlo 
algorithm and a powerful finite-size-scaling extrapolation method. For q = 3 we 
obtain good control up to correlation length £ ~ 5000; the data are consistent 
with = Ae 2l3 (3 p (l + a\e~^ + . . .) as f3 — > oo, with psl, The staggered 
susceptibility behaves as Xstagg ~ £ 5 ^ 3 . For q = 4 the model is disordered (£ < 
2) even at zero temperature. In appendices we prove a correlation inequality 
for Potts antiferromagnets on a bipartite lattice, and we prove ergodicity of the 
WSK algorithm at zero temperature for Potts antiferromagnets on a bipartite 
lattice. 
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1 Introduction 



The Potts model []l], |^, |3j plays an important role in the general theory of critical 
phenomena, especially in two dimensions ||, ^ ||, and has applications to various 
condensed- matter systems 0. Ferromagnetic Potts models have been extensively 
studied over the last two decades, and much is known about their phase diagrams 
0. and critical exponents |5], |6|, 0]. But for antiferromagnetic Potts models, many 
basic questions remain open: Is there a phase transition at finite temperature, and 
if so, of what order? What is the nature of the low-temperature phase(s)? If there 
is a critical point, what are the critical exponents and the universality classes? Can 
these exponents be understood (for two-dimensional models) in terms of conformal 
field theory? 

One thing is known rigorously || |J: for q large enough (how large depends on the 
lattice in question), the antiferromagnetic g-state Potts model has a unique infinite- 
volume Gibbs measure and exponential decay of correlations at all temperatures, 
including zero temperature: the system is disordered as a result of the large ground- 
state entropy. However, for smaller values of q, phase transitions can and do occur. 
Moreover, for these antiferromagnetic models the nature of the phase transition is 
highly lattice-dependent, in sharp contrast to the universality typically enjoyed by 
ferromagnets. Thus, one expects that for each lattice £ there will be a value q c {£) 
such that 

(a) For q > q c (C) the model has exponential decay of correlations uniformly at all 
temperatures, including zero temperature. 

(b) For q = q c (C) the model has a critical point at zero temperature. 

(c) For q < q c (C) any behavior is possible. Often (though not always) the model 
has a phase transition at nonzero temperature, which may be of either first or 
second order. 

The problem, for each lattice, is to find q c (C) and to determine the precise behavior 
for each q < q c {C). 

In this paper we report the results of a large-scale Monte Carlo study of the 3-state 
and 4-state antiferromagnetic Potts models on the (two-dimensional) square lattice, 
using the Wang-Swendsen-Kotecky (WSK) [II], O, [T2j cluster algorithm.]] We use a 
powerful finite-size-scaling (FSS) extrapolation method |14[ [15], [U| [17|, [0| [19], |(], [21], 



|22|j to estimate the infinite-volume correlation length £ and staggered susceptibility 
Xstagg- Using lattices up to 1536 x 1536, we can attain an accuracy of a few percent 
on £ and Xstagg at correlation lengths £ as large as 5000. This allows us to conjecture 
the exact form of the critical behavior for the 3-state model. 

The g-state Potts model is defined by the reduced Hamiltonian 

« = ~Jj2 6 ^y > (1-1) 



1 A preliminary version of this work has appeared previously [[l3[ 
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where the sum runs over all nearest-neighbor pairs of lattice sites, and each spin takes 
values a x G {1,2, ... ,q}. The antiferromagnetic case corresponds to J = —[3 < 0. 
Henceforth we restrict attention to the model ( |1 . 1|) on the square lattice. 

Baxter [B, ^] has determined the exact free energy (among other quantities) for 
the square-lattice Potts model on two special curves in the (J, g)-plane: 

e J = l±y/q (1.2) 
e J = -l±J±-q (1.3) 



Curve ( |l-4 f ) 1S known to correspond to the ferromagnetic critical point, and Baxter 
p3| conjectured that curve ([L3[+) corresponds to the antiferromagnetic critical point. 
For q = 2 this gives the known exact value |24| ; for q = 3 it predicts a zero-temperature 
critical point (J c = — oo), in accordance with previous belief [^ |26]Q; and for q > 3 
it predicts that the putative critical point lies in the unphysical region (e Jc < 0), so 
that the entire physical region — oo < J < lies in the disordered phase. In other 
words, Baxter p3| predicts that q c (square lattice) = 3, a prediction that we will verify 
numerically in this paper. 

Some properties of the zero-temperature critical point for q = 3 are known (non- 
rigorously) as a consequence of the mapping of this model onto a height model, whose 



long-wavelength behavior is that of a massless Gaussian [30, 31, 32. 33. 34]. In partic- 
ular, the critical exponents associated to the staggered and uniform magnetizations 
are predicted |3C], [32] to be r) stagg = 1/3 and r\ u = 4/3, respectively.^ However, the 
approach to the critical point is much less well understood. In what way, for example, 
do £ and Xstagg diverge as (3 — > oo? 



Nightingale and Schick [35], using a phenomenological-renormalization method 
based on infinite strips of width 2-8, claimed that the correlation length diverges as 
£ ~ exp(c/3~ L3 ). Wang, Swendsen and Kotecky JK| [H| , using Monte Carlo, claimed 



to confirm this latter behavior. But this behavior seems a priori implausible to us: 
the fundamental variable in the Potts model is t = e J , so an ordinary power-law 
critical point £ ~ (t — t c )~ u with t c = would correspond to £ ~ e u/3 . Moreover, we 
suspect that this model can be exactly solved (at least in the sense of determining 
the exact asymptotic behavior as (3 — > oo), in which case v would most likely be a 
rational number. We are unable to imagine any mechanism leading to £ ~ exp(c(3 K ) 
with 1. 

In this paper we shall present numerical evidence that strongly suggests the asymp- 
totic behavior 

£(P) = Ae w f3 p [l + ai e~ p + a 2 e~ w + . . .] (1.4) 

with p m 1. The critical exponent v = 2 found here corresponds to an operator with 
scaling dimension X — 2 — l/i/ = 3/2, which is one of the possibilities proposed by 

2 Note also that the q — 3 model is exactly soluble at zero temperature in an arbitrary magnetic 
field |26|, ^9) ; this might increase one's suspicions that the zero-temperature zero- field case is 
critical. 

3 These predictions have recently been verified numerically to high precision [^3 . 
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Saleur |36], p. 248] — though not the one he considered most likely! The multiplicative 
logarithmic correction (3 P ~ {\ogt) p is harder to understand theoretically. Indeed, 
our data can arguably be reconciled with p = at the price of including additive 
corrections to scaling based on fractional powers of e _/3 : 

£(/3) = Ae 2/3 [l + a ie ~ Al/3 + a 2 e~ A2/3 + . . .] (1.5) 

with Ai ~ 0.5. But this Ansatz too has its theoretical difficulties; see Section |7.1| 
for discussion of all these issues. We hope, in any case, that the numerical results 
presented here will serve as useful clues toward the exact solution of this model. 



The works of Saleur |37], |36| and Henley |38| provide some tantalizing ideas in this 
direction. 

As for the staggered susceptibility, the predicted critical exponent r\ = 1/3 |30], |32| 
leads via the scaling law 7 = (2 — r\)v to the behavior Xstagg ~ £ 5//3 - Our numerical 
results are consistent with this prediction, unmodified by any further powers of (3. 

The plan of this paper is as follows: In Section [| we set the notation and recall 
briefly our finite-size-scaling extrapolation method and the Wang-Swendsen-Kotecky 
(WSK) Monte Carlo algorithm. In Section |3] we report our raw data. In Section |] we 
analyze our static data for the 3-state model, using the finite-size-scaling extrapolation 
method. In Section |5] we analyze the dynamic critical behavior of the WSK algorithm 
for the 3-state model. In Section |6] we analyze the data for the 4-state model. In 
Section [7| we summarize our conclusions and discuss prospects for future work. In 
Appendix |A| we prove a correlation inequality for antiferromagnetic Potts models on 
a bipartite lattice. In Appendix [TJ we prove the ergodicity of the WSK algorithm at 
zero temperature for antiferromagnetic Potts models on a bipartite lattice. 



2 Preliminaries 

2.1 Definitions and Notation 

The g-state Potts model is defined by the reduced Hamiltonian 

where the sum runs over all nearest-neighbor pairs of lattice sites, and each spin takes 
values a x G {1, 2, . . . , q}. The antiferromagnetic case corresponds to J = —(3 < 0. It 
is useful to represent the q possible values of the spin a x by unit vectors e^, . . . , G 
1R 9_1 pointing from the center to the vertices of a hypertetrahedron; these vectors 
satisfy 

eW • eW = = { x :e , , (2.2) 

We denote this "vectorial" spin by cr x = e^ ax \ 
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The two-point correlation function G(x, y) is defined by 



G(x,y) = <<r, ■ tr v ) = - - \ . (2.3) 

Henceforth we exploit translation invariance and write G(x,y) = G(x — y). We also 
define the Fourier-transformed correlation function at wavevector ("momentum") p: 

G(p) = J2e ip ' x G(x). (2.4) 

X 

On the square lattice, the relevant staggering wavevector for antiferromagnetic Potts 
models is 

Pstagg = (tT, 7r) , (2.5) 

in the sense that G(p) is maximum at p — (ir,7r): this follows from the correlation 
inequality proven in Appendix [A|. On a finite L x L lattice with periodic boundary 
conditions, we also define the four smallest nonzero wavevectors, 

Pmin.il = (±2tt/L,0) (2.6a) 
p min , ±2 = (0,±2tt/L) (2.6b) 

We wish to study the following quantities: 

(a) The energy^ 

E = (<We> , (2-7) 

where e stands for any nearest neighbor of the origin. 

(b) The staggered susceptibility 

Xstagg = G(p sta gg) ■ (2.8) 

Note that on a finite lattice this is well-defined only if the lattice size L is even. 

(c) The second-moment correlation length, which is defined in finite volume by 

c [(Xstagg/ F s tagg) ~ 1] ^ /„ „x 

4l = . / /r x > (2.yj 

2 sm(7r/L) 

where 

Fstagg = G(p s tagg + Pmin,±l) = G(p s tagg + Pmin,±2) (2-10) 

is the correlation function at the wavevectors closest to Again L must be even. 

4 Here E is the mean energy per link in the antiferromagnetic model; we have chosen this 
normalization in order to have < E < 1, with E = for an antiferromagnetic ground state and 
E = 1 for a ferromagnetic ground state. The mean energy per site is of course 2E. 

5 See e.g. |3!| equations (4.11)-(4.13)] for the definition of the second-moment correlation length 
in a ferromagnetic model, along with its motivation. Here we make the obvious transcription to an 
antiferromagnetic model that orders at momentum Pstagg — (tTj tt). 
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All these quantities can be expressed as expectations involving the following ob- 
servables: 



stagg 



3~stagg 

£ = 
Thus, we have 



^] e lp stagg x (Tx j 

^ y fD^(Pstagg~^Pmin,+l)' x ( 



+ 



e i (Pstagg+Pmin,+2 



(2.11a) 

(2.11b) 
(2.11c) 



Xstagg 
Fstagg 

E 



1 

V 



i?{M 2 stagg ) 



y {3~stagg) 



(2.12a) 
(2.12b) 
(2.12c) 



where V = L 2 is the number of lattice sites. 

In addition to studying the (static) behavior of the antiferromagnetic Potts model, 
we are also interested in studying the dynamic critical behavior of the Wang-Swendsen- 
Kotecky ( WSK) Monte Carlo algorithm. So let A be an observable (i.e. a function of 
the spin configuration {a}). We define the unnormalized autocorrelation function 



C AA (t) = (A s A s+t ) - (A) 2 , 



(2.13) 



where expectations are taken in equilibrium, and the corresponding normalized au- 
tocorrelation function 

PaaH) = C AA (t)/C AA (0) . (2.14) 
We furthermore define the integrated autocorrelation time 



oo 

T~int,A = ~ XI PAA(t) 
Z t=-oo 



1 



Y,PAA(t) ■ 



(2.15a) 
(2.15b) 



t=i 



[The factor of \ is purely a matter of convention; it is inserted so that Ti nttA ~ r if 
PAA(t) ~ e _ l'l/ r with r >> 1.] Finally, we define the exponential autocorrelation times 



and 



Texp A = lim SUp j— — - 

t^oo -log \p AA [t)\ 



l~exp SUp T eX p jA . 

A 



(2.16) 
(2.17) 
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Note that r exp = r exP:A whenever the observable A is not orthogonal to the slowest 
mode of the system. 

The integrated autocorrelation time controls the statistical error in Monte Carlo 
measurements of (A). More precisely, the sample mean 

in 

A = (2.18) 



has variance 



1 n 

var(A) = - E C AA (r-s) (2.19a) 



rr 



r,s = 1 



1 ^ It 



E (!- U )C^(t) (2.19b) 
(2r intA ) C AA {0) for n > r (2.19c) 



i = — (n— 1) 

1 

n 



Thus, the variance of A is a factor 2rj nt>J 4 larger than it would be if the {^4f} were 
statistically independent. Stated differently, the number of "effectively independent 
samples" in a run of length n is roughly n/2r inttA . The autocorrelation time r inttA 
(for interesting observables A) is therefore a "figure of (de)merit" of a Monte Carlo 
algorithm. 

The integrated autocorrelation time Ti n t,A can be estimated by standard proce- 
dures of statistical time-series analysis |^0|, These procedures also give statisti- 



cally valid error bars on (A) and Ti nt>A . For more details, see f42| , Appendix C] or [43 
Section 3]. In this paper we have used a self-consistent truncation window of width 
CTint,A, where c = 6; this choice is reasonable whenever the autocorrelation function 
PAAif) decays roughly exponentially, as it does here (see Section ^]2| below). 

In setting the error bars on £ [defined in ( p.9|) 1 we have used the triangle inequality; 
such error bars are overly conservative, but we were too lazy to measure the cross- 
correlations between M. 2 stagg and T s tagg- (This was a mistake, and in future work we 
will measure these cross-correlations.) 



2.2 Finite-Size-Scaling Extrapolation Method 

In the theory of critical phenomena we are usually interested in infinite-volume 
systems, but Monte Carlo simulations are perforce carried out on finite lattices. One 
traditional approach has been to run on lattice sizes L > 6£, which are large enough so 
that the finite-size corrections are negligible. In the past few years, however, methods 
have become available for extrapolating finite-size data to L = oo, based on finite-size- 
scaling (FSS) theory; these methods allow one to work, for a given lattice size L, at 
correlation lengths £ much larger than were previously attainable. In this subsection 
we review an extremely powerful and general method of this kind, due originally to 
Luscher, Weisz and Wolff |L4| and more recently elaborated by our group ]20], [21], p2| . 
This extrapolation method plays a crucial role in the present work, as it allows us 
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to reach correlation lengths £ of order 5000, whereas in the traditional approach we 
would have been limited to £ < 250. 

Consider, for simplicity, a model controlled by a renormalization-group fixed point 
having one relevant operator. (The parameter (3 may be multidimensional, but only 
one direction in /9-space should be relevant in the RG sense. In other words, the con- 
tinuum limit should be unique modulo length rescalings.) Let us work on a periodic 
lattice of linear size L. Let L) be a suitably defined finite- volume correlation 
length, such as the second-moment correlation length defined by ( |2.9| ); and let O be 
any long-distance observable (e.g. the correlation length or the susceptibility). Then 
finite-size-scaling theory [[0], [|5|, predicts that 



0(P,oo) 



/ (£(/?, oo)/L) + 0(C^L^) , (2.20) 



where fo is a universal (though usually unknown) function and u is a correction-to- 
scaling exponent.[] It follows that if s is any fixed scale factor (usually we take s = 2), 
then 

= F (t(/3,L)/L; s) + 0(C",L-) , (2.21) 

where Fq can easily be expressed in terms of fo and fa. (Henceforth we shall suppress 
the argument s if it is clear from the context.) In other words, if we make a plot of 
0(P, sL)/0(/3, L) versus £(/?, L)/L, then all the points should lie on a single curve, 
modulo corrections of order and L~ w . 

Our method proceeds as follows^: Make Monte Carlo runs at numerous pairs (j3, L) 
and (/3,sL). Plot (D(j3, sL)/0(j3, L) versus L)/L, using those points satisfying 
both L) > some value £ m j n and L > some value L m i n . If all these points fall 
with good accuracy on a single curve — thus verifying the Ansatz Q2.21|) for £ > £ m j n , 

6 This form of finite-size scaling assumes hyperscaling, and thus is expected to hold only below 
the upper critical dimension of the model. See e.g. Jl6|, Chapter I, section 2.7]. Note also that when 
we say fo is "universal" , we mean only that it is the same for all models in a given universality 
class. Of course fo varies from one universality class to another. 

7 Our method ||(J |l], ||| is essentially identical to that of Liischer, Weisz and Wolff @. The 
principal difference is that Liischer et al. choose carefully their runs (/?, L) so as to produce only a few 
distinct values of x = £(/3, L)/L, while we attempt to cover an entire interval of x. Which approach 
is preferable depends on one's aims and on the available CPU time. Also, the motivations are 
somewhat different: the primary aim of Liischer et al. [ fTi) is to compare the asymptotic behavior of 
the finite-size-scaling functions Fo(x) at large x to the perturbative predictions; they did not discuss 
the possibility of obtaining extrapolations to L = oo at each fixed (3, although this is of course 
implicit in their method. The method of Kim |15, [l6[ [l7| [f8[ |l9) is also very closely related, but he 



compares lattice size L to oo rather than to sL; this is a (slight) disadvantage. Our method also 
has many features in common with that used by Flyvbjerg and Larsen fl7L [f8| to extrapolate their 
1/iV-expansion finite-lattice data. It should be emphasized that all these methods are completely 



general; although they were historically first applied to asymptotically free theories 1 14 , they are in 
no way limited to this case. Note also that all these methods share the property of working only 
with observable quantities (£, O and L) and not with bare quantities (/?). Therefore, they rely only 
on "scaling" and not on "asymptotic scaling" ; and they differ from other FSS-based methods such 
as phenomenological renormalization p9[. 
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L > L m i n — choose a smooth fitting function Fq. Then, using the functions F% and 
Fqi extrapolate the pair (£, O) successively from L — > sL — > s 2 L oo. 
We have chosen to use functions -F^ of the formF] 

F (x) = 1 + a ie ~ 1/x + a 2 e- 2/:c + . . . + a n e- n/x . (2.22) 

(Other forms of fitting functions can be used instead.) This form is partially motivated 
by theory, which tells us that in some cases Fq(x) — > 1 exponentially fast as x — > 0. 
Typically a fit of order 3 < n < 13 is sufficient; the required order depends on the 
range of x values covered by the data and on the shape of the curve. Empirically, we 
increase n until the x 2 of the fit becomes essentially constant. The resulting \ 2 value 
provides a check on the systematic errors arising from corrections to scaling and/or 
from the inadequacies of the form ( |2.22j ) . 

The statistical error on the extrapolated value of C 00 (/5) = 0(/3, oo) comes from 
three sources: 

(i) Error on 0(j3, L), which gets multiplicatively propagated to O^. 

(ii) Error on L), which affects the argument x = £(ft, L)/L of the scaling func- 
tions F% and Fq. 

(iii) Statistical error in our estimate of the coefficients ax, . . . , a n in F% and Fq. 

The errors of type (i) and (ii) depend on the statistics available at the single point 
(ft, L), while the error of type (iii) depends on the statistics in the whole set of runs. 
Errors (i) + (ii) [resp. (i) + (ii) + (iii)] can be quantified by performing an auxiliary Monte 
Carlo experiment in which the input data at (ft, L) [resp. the whole set of input data] 
are varied randomly within their error bars and then extrapolated (we call this the 
method of "fake data sets").[] 

8 In performing this fit, one may use any basis one pleases in the space spanned by the functions 
{e~ k / x }\<k<n] the final result (in exact arithmetic) is of course the same. However, in finite-precision 
arithmetic the calculation may become numerically unstable if the condition number of the least- 
squares matrix gets too large. In particular, this disaster occurs if we use as a basis the monomials 
t k (where t — e~ 1 l x ). The trouble is that these monomials are "almost collinear" in the relevant 
Hilbert space L 2 (n) defined by fi(t) = Wi S(t — ti), where U are the values of t = e^ L ^^- L ^ arising 
in the data pairs and Wi — 1/ [error on 0(2L)/0(L)) 2 are the corresponding weights. To avoid this 
disaster, we should seek to use a basis that is closer to orthogonal in i 2 (/i). Of course, exactly 
orthogonalizing in L 2 (/i) is equivalent to diagonalizing the least-squares matrix, which is unfeasible; 
but we can do well enough by using polynomials with zero constant term that are orthogonal with 
respect to the simple measure w(t) = t a (t max — t) b on [0, t max ], where a and b are some chosen 
numbers > — 1. These polynomials are Jacobi polynomials fk(t) = tPj: b ^ L +2 \2t/t max — 1) for 
1 < k < n pp. 321-328]. The idea here is that the measure w(t) = t a (t max — t) b should roughly 
approximate the measure fJ-(t). We have here used a = 0, b = —3/4; but the performance is very 
insensitive to the choices of a and b. This cleverness in the choice of basis vastly improves the 
numerical stability of the result, by reducing the condition number of the matrix arising in the fit. 
Typical condition numbers using Jacobi polynomials are ~ 20 for n — 3 and ^100 for n = 9. 

9 In principle, £ and O should be generated from a joint Gaussian with the correct covariance. 
We ignored this subtlety and simply generated independent fluctuations on £ and O. 
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The discrepancies between the extrapolated values from different lattice sizes at 
the same j3 — to the extent that these exceed the estimated statistical errors — can 
serve as a rough estimate of the remaining systematic errors. More precisely, let Oi 
(i = 1, . . . ,m) be the extrapolated values at some given j3, and let C = (C^)^- =1 
be the estimated covariance matrix for their statistical errors.0 [Errors of type (iii) 
induce off-diagonal terms in C] Then we form the weighted average 

O = (t {C-%o\ I ( ± (C-%) , (2.23) 
V<i=i / / \m=i / 

the error bar on the weighted average 

/ m X- 1 / 2 

* = iE(C-%\ , (2.24) 
and the residual sum-of-squares 

m 

K = E (Oi - 0){C-%{O 3 - O) . (2.25) 

Under the assumptions that 

(a) the fluctuations among the Oi, . . . , O m are purely statistical [i.e. there are no 
systematic errors in the extrapolation], and 

(b) the statistical error bars are correct, 

TZ should be distributed as a \ 2 random variable with m — 1 degrees of freedom. 
Moreover, the sum of 1Z over all the values of j3 should be distributed as a x 2 random 
variable with — degrees of freedom.^] In this way, we can search for values of 
(3 for which the extrapolations from different lattice sizes are mutually inconsistent; 
and we can test the overall self-consistency of the extrapolations. 

A figure of (de)merit of the method is the relative variance on the extrapolated 
value C?oo(/3), multiplied by the computer time needed to obtain it.0 We expect this 
relative variance-time product [for errors (i) + (ii) only] to scale as 

RVTP(/?,L) « U{P) d+Zint '°G {U{P)/L) , (2.26) 



10 This covariance matrix is computed from the auxiliary Monte Carlo experiment mentioned in 
the preceding paragraph. Since this C is only a statistical estimate, the values of O, a and 1Z will 
vary slightly from one analysis run to the next. 

11 This latter statement is not quite correct, as it ignores the correlations between the various Oi 
at different (3, which are induced by errors of type (iii). [Correlations between different Oi at the 
same (3, which are also induced by errors of type (iii), are included in (|^23|)-(^25|).1 

12 At fixed ((3, L), this variance-time product tends to a constant as the CPU time tends to infinity. 
However, if the CPU time used is too small, then the variance-time product can be significantly larger 
than its asymptotic value, due to nonlinear cross terms between error sources (i) and (ii). 
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where d is the spatial dimension and Zi n t,o ls the dynamic critical exponent of the 
Monte Carlo algorithm being used; here Go is a combination of several static and 
dynamic finite-size-scaling functions, and depends both on the observable O and on 
the algorithm but not on the scale factor s. As £,oo/L tends to zero, we expect Go to 
diverge as (£oo/ ' L)~ d (since it is wasteful to use a lattice L ^> £00). As £oo/£ tends to 
infinity, we expect Go ~ (£,oo/L) p for some power p (see |22[] for details). Note that 
the power p can be either positive or negative. If p > 0, there is an optimum value 
of £00/ L; this determines the best lattice size at which to perform runs for a given f3. 
If p < 0, it is most efficient to use the smallest lattice size for which the corrections 
to scaling are negligible compared to the statistical errors. [Of course, this analysis 
neglects errors of type (iii). The optimization becomes much more complicated if 
errors of type (iii) are included, as it is then necessary to optimize the set of runs as 
a whole.] 

The reader is referred to [[H], |22] for a fuller treatment of this extrapolation method, 



in particular the finite-size-scaling theory and the analysis of the propagation of sta- 
tistical errors. 

Let us make one final comment about the physics contained in the scaling function 
F^(x). At the critical point /3 C , the correlation length £(/3 c , L) is proportional to L: 
one thus has 

lim — - = some value x* (2.27) 

and 

Therefore, x* is determined by the relation 

F^(x*;s) = s . (2.29) 

The constant x* is characteristic of the massless field theory corresponding to the 
given critical point, on a continuum torus with aspect ratio 1; for two-dimensional 
models it should in principle be calculable via conformal field theory. Likewise, for 
any observable O which behaves in the critical region like O ~ one has 

F (x*; s) = s lo/u . (2.30) 



2.3 Wang-Swendsen-Kotecky (WSK) Algorithm 

About a decade ago, Wang, Swendsen and Kotecky (WSK) |IJ [TT[] proposed an 
elegant and extraordinarily efficient Monte Carlo algorithm for simulating the anti- 
ferromagnetic g-state Potts model on an arbitrary finite graph G. The elementary 
update of their algorithm goes as follows: Choose at random two distinct "colors" 
i,j G {1, . . . , q}; freeze all the spins a x currently taking values k 7^ and allow the 
remaining spins to take value either i or j. The induced model is thus an antiferro- 
magnetic Ising model (in zero magnetic fieldQ) on a subgraph of G; this model can 

13 The key point here is that the interaction energies i—k and j—k are equal. This guarantees that 
the induced Ising model has zero magnetic field. 
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be updated by any legitimate Monte Carlo algorithm, such as the Swendsen-Wang 
algorithm [[51] or Wolff's single-cluster variant |52| (we use the former). It is easy to 
see that this algorithm leaves invariant the Gibbs measure of the underlying Potts 
model. 

At zero temperature (/3 = oo) the antiferromagnetic Potts model reduces to the 
equal-weight distribution on g-colorings of G, and the WSK algorithm becomes: in- 
dependently for each connected cluster of i—j spins, either leave that cluster as is or 
else flip it (interchanging i and j). 

The WSK algorithm is trivially seen to be ergodic at any nonzero temperature. 
However, the ergodicity at zero temperature is a very subtle problem, which has thus 
far been only partially resolved. Lubin and Sokal |TI| showed that for q = 3 the 



WSK algorithm is non- ergodic at zero temperature on periodic square lattices of size 
3m x 3n where m and n are relatively prime. On the other hand, in Appendix [B] 
we shall prove that the WSK algorithm is ergodic for all q whenever the graph G is 
bipartite; in particular, this happens on periodic square lattices of size mxn whenever 
m and n are both even. For other cases the ergodicity is an open problem. 

It should be noted that the problem of ergodicity at zero temperature is not merely 
a theoretical one, even if all our simulations are performed at nonzero temperature. If 
the algorithm is non-ergodic at (3 = oo, then the autocorrelation time must diverge as 
(3 — > oo, even on a fixed finite lattice.™ This behavior, if it occurs, could be a severe 
impediment to high-precision Monte Carlo study of the model. So it is fortunate that 
our model here lies precisely in the situation for which ergodicity has been proven: 
the periodic square lattice with L even. 



3 Summary of our Runs 

We simulated the square-lattice Potts antiferromagnets for q = 3 and q — 4, using 
the WSK algorithm with standard (multi-cluster) Swendsen-Wang updates of the 
induced Ising model. 

For q = 3 we ran on LxL periodic lattices with L = 32, 64, 128, 256, 512, 1024, 1536 
at 149 different pairs (/5, L) in the range 2.0 < (3 < 6.0 (corresponding to 5 ^ £oo ~ 
20000). Our data cover rather densely the range 0.09 < x = C,(/3, L)/ L < 0.63. Each 
data point comprises between 2 x 10 5 and 2.2 x 10 7 iterations of the WSK algorithm, 
which corresponds to anywhere from 40000r to 5 x 10 6 r. We discarded the first 
10000 iterations of each run, which ought to be more than enough for equilibration 
(> 2000t!); we also made spot checks for evidence of initialization bias after the 
discard interval, and found none. Most of the runs used a random initial configu- 
ration ("hot start"). In some cases with L = 32,64 we made multiple independent 

14 More precisely, it is the exponential autocorrelation time T exp (corresponding to the slowest 
mode in the system) which must diverge. The integrated autocorrelation time Ti n t,A for any given 
observable A need not diverge; that depends on the choice of A. Moreover, even if it does diverge, 
the divergence could be very weak, if A has "weak overlap" with the slowest mode. On the other 
hand, the divergence of T exp already calls into question the convergence to equilibrium, by raising 
the specter of "metastability" . 
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runs, some of which used an antiferromagnetically ordered initial configuration ( "cold 
start"); statistical tests showed complete agreement between the runs. The data from 
independent runs were merged statistically in the usual way. The raw data for q = 3 
can be found in Table [I]. 

The CPU time for our program is 4.4L 2 /zs/sweep on an IBM RS-6000/370, and the 
total CPU time was about 2.5 years on this same machine. (This is an "equivalent" 
figure: in fact our runs were performed on a variety of mostly slower machines in both 
New York and Belo Horizonte, so the actual total CPU time was more than this.) 

For q = 4 it sufficed to make a small number of runs for L = 32, 64; the total CPU 
time was less than 3 days on an IBM RS-6000/370. The raw data can be found in 
Table §. 



4 Data Analysis: q = 3, Static Quantities 

In this section we analyze the static data for q = 3. Concerning the correlation 
length £ and the staggered susceptibility Xstagg, we first extrapolate these quantities 
to infinite volume (Section |4.1|) and then analyze the behavior as (3 — > oo of the 
extrapolated data (Section |4~2{ ). We conclude by taking a brief look at the energy E 
(Section Q. 



4.1 Extrapolation to Infinite Volume 

We shall extrapolate £ and Xstagg to infinite volume using the method of Section 
|2.2| with scale factor s = 2. This method is specified by three parameters: the cut 
points ^min and L m i n , and the interpolation order n. We shall therefore vary these 
parameters systematically and study the systematic errors attributable to them. 



4.1.1 Correlation Length 

In Table |3] we report the quality of the fit — chi-squared (x 2 )> number of degrees 
of freedom (DF), x 2 /DF, and the corresponding confidence level[^] — for the function 
F$ = £(/3, 2L) /£(/?, L) as a function of the interpolation order n and the cut point 
L min ] here we have used £ m j n = 10. (We tried also £ m j n = 20 and the results are 
virtually unchanged.) A reasonable x 2 is obtained when n > 5 and L m i n > 64; 
a slightly better x 2 /DF is obtained by taking L min = 128. Further increases in 
n and/or L min do not improve the x 2 /DF.[^] Our preferred fit is therefore n = 5, 

15 "Confidence level" is the probability that \ 2 would exceed the observed value, assuming that 
the underlying statistical model is correct. An unusually low confidence level (e.g. less than 5%) thus 
suggests that the underlying statistical model is incorrect. Here this may be due to an inadequate 
interpolation Ansatz (too low n) or to corrections to scaling (too low £ m i„ or L m i n ). Another possible 



cause of unusually low confidence levels will be discussed in Section 4.1.2. 

16 Indeed, for L m i n = 256 the x 2 /DF is worse than for L m i n — 64, 128; it is, in fact, as bad as for 
Lmin — 32! We do not understand the reason for this behavior, which may be simply a statistical 
fluctuation. 
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imin = 10 and L min = 128: we get 

F ( (x) = 1 + 0.896737e~ 1/x + 29.243141e~ 2/:c - 253.811947e~ 3/:!; 

+ 2092.996892e" 4/:r - 5334.794958e~ 5/a; , (4.1) 

which is plotted in Figure |I[ This fit is reliable only in the interval where there are 
data points contributing to it, namely x < x max ~ 0.629781. 

We remind the reader that our raw-data error bars on £ are overestimates, as 
a result of our use of the triangle inequality; therefore, the x 2 values reported in 
Table |3| are underestimates , and only their relative magnitudes can be considered to 
be reliable. The absolute quality of the fit is, therefore, not as good as it looks. This 
will be further discussed below in connection with Xstagg 

(Section gXp . 

We can make a crude estimate of the universal value x* defined by F^(x*) = 2 [cf. 
O)]. On the one hand, F^(x max ) = 1.987350, so x* > x max « 0.629781. On the 



other hand, F^(x max ) = 3.080022, so if we extrapolate linearly for x > x max , we get 
x* = 0.633888. This is not far from the value x* = 0.633983 obtained by taking ( |4.1| ) 
seriously even for x > x max . So it is a fair guess that x* ~ 0.633888.Q 

In Figure [|(a) we plot the deviations from our preferred fit together with their 
error bars. The points with L = 32 show weak (< 0.01) but apparently statistically 
significant deviations, of positive sign, in the interval 0.52 < x < 0.625: see the 
blow-up of this region in Figure @(b).0 It is a reasonable guess that these deviations 
arise from systematic corrections to scaling. However, a careful test of this hypothesis 
would require higher statistics and more densely spaced points; in addition, it would 
be useful to obtain data with extremely high statistics (~ 10 8 sweeps) on smaller 
lattices than we have bothered to use here (L = 16 and even L = 8), in order to 
observe a stronger correction-to-scaling "signal". The points with L > 64 do not 
seem to show any systematic deviations from the fit. Because we observe statistically 
significant corrections to scaling on only one lattice size, we are unable to make any 
firm statement about the L-dependence of the correction-to-scaling term (which we 
expect to be of the form L~ A , where A > is a correction-to-scaling exponent). All 
we can say is that A is not too small, since otherwise the correction to scaling would 
be observable also on the L = 64 lattice. But we are unable to say whether, for 



example, A fa 1 or A rs 2. Let us remark that Salas and Sokal [31| have recently 
predicted, on the basis of the height representation of the zero-temperature model 
|50], |3T], ^ [53] g, that A = 2 (at least when (5 = oo). 

We now compute the extrapolated values ^ , using the function F^ given by our 
preferred fit as well as by several alternative fits that use more or less stringent choices 
of L m i n . (In all cases we take n = 5; the results for n = 6 or n = 7 differ in almost all 
cases by < 0.4 standard deviations.) The statistical error bars on ^ are computed 
by an auxiliary Monte Carlo process using "fake data sets", as described in Section 



17 After completion of this work, Salas and Sokal j34| obtained x* = 0.63457 ± 0.00033 by high- 
precision simulation of this model at (3 — oo. 

18 The L = 32 point at x 0.633 lies outside the range of x covered by the fit (x max ~ 0.630 
128), so the negative deviation exhibited by this point may not be meaningful. 
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2~2l , and include errors of types (i) + (ii) + (iii). In particular, the correlations between 
the extrapolated values ^ from different lattice sizes at the same (3 (though not 
those at different (3 values) are taken account of in this computation. For each (3 
we compute the weighted average ( |2.23| ) of the various estimates £oo, along with the 
statistical error bar (|2.24|) and the 1Z value ( |2.25| ). These are reported in Table f|. 
The statistical errors on are of order 1% (resp. 2%, 3%, 5%) at £oo ~ 1000 (resp. 
2000, 5000, 10000). These extrapolated values from different lattice sizes at the same 
f3 are found to be consistent within statistical errors: only two of the 42 (3 values has 
an 1Z value too large at the 5% level; and summing all [3 values we have ^TZ — 43.45 
(75 DF, level = 99.9%). This unusually low chi-squared is probably due, at least 
in part, to our overestimation of the raw-data error bars on which leads to an 
overestimation of the error bars on l;^. 

The discrepancies between the extrapolations with L min = 64, 128, 256 are in 
almost all cases less than half the quoted statistical error.0 We are thus reasonably 
confident that we have obtained quantitative control over the systematic errors due 
to corrections to scaling, and that their effect can be at most to double the quoted 
statistical errors. 



4.1.2 Staggered Susceptibility 

Next we carry out an analogous analysis for the staggered susceptibility Xst agg - 
Note that for this observable our raw-data error bars are reliable, so the absolute x 2 
values can be taken seriously. 

In Table [5] we report the quality of the fit for the function F Xstagg = Xsta gg ({3, 21/) / Xst agg (P, L) 
as a function of the interpolation order n and the cut point L m i n ; here we have used 
imin = 10. (We tried also £ min = 20 and the results are virtually unchanged.) The 
% 2 /DF is smallest when n > 6 and L min > 128.Q However, this x 2 /DF is ~ 2, rather 
than the ~ 1 that it ought to be; and as a result, the confidence levels are extremely 
low (of order 10~ 6 ). We do not understand this behavior, but we can make a few 
observations: 

(a) Clearly, the explanation cannot be either an inadequate fitting function or 
corrections to scaling, because increasing n and/or L min and/or £ min does not improve 
the fit. 

(b) One possible explanation might be that our raw-data error bars are underesti- 
mated; but we tried various alternative statistical methods, such as breaking up long 
runs into sub-runs, and all gave compatible error bars. So we do not think that this 
is the problem. 

(c) It is worth noting that nearly all the points with poor x 2 come from the region 
of large x, particularly x ^ 0.59. Indeed, if we look separately at the contributions 

19 On the other hand, the extrapolations with L m j n = 32 frequently deviate from our preferred 
Lmin — 128 extrapolation by as much as ~ la. These deviations are probably a correction-to-scaling 
effect on the borderline of statistical significance. 

20 The x 2 /DF gets worse for L m i n = 256, just as it does for the correlation length (see footnote 
|l6| above) . 
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to x 2 (using the fit with n — 6, £ m j„ = 10, L min = 128) coming from the intervals 
x < 0.53, 0.53 < x < 0.59 and x > 0.59, we find that the first interval contains 29 
points which together contribute 24.31 to the x 2 (x 2 /DF = 0.84, level = 71.4%), the 
second interval contains 10 points which together contribute 16.65 to the x 2 (x 2 /DF 
= 1.67, level = 8.2%), while the third interval contains 32 points which together 
contribute 90.49 to the x 2 (x7 DF = 2.83, level = 1.7 x HT 7 ). 

(d) This behavior can also be seen in Figure |3], where we show the deviations from 
the fit n — 6, £ m j n = 10, L min = 128 together with their error bars. The points with 
L = 32 show weak (< 0.02) but statistically significant deviations, of positive sign, 
in the interval 0.53 ^ x ^ 0.625: see the blow-up of this region in Figure |3|(b).0 
However, for L > 64 the corrections to scaling have become completely invisible. On 
the other hand, for x > 0.58 [see the enlarged view in Figure |3|(c)] we see that several 
points deviate from the fitting curve by 2-4 standard deviations; but there does not 
seem to be any systematic trend to these deviations, nor do the absolute deviations 
seem to be larger for smaller L. So, once again, it is unlikely that these deviations are 
caused by corrections to scaling, or indeed by any process that causes a systematic 
bias. 

(e) One possible cause of the unusually high x 2 is the following: We treated 
the points in the fit - - which of course correspond to pairs ((3, L) / ((3, 2L) - - as 
statistically independent; but this is not quite right, as the same raw-data point 

2L) can contribute to two pairs, namely (/3, L)/(j3, 2L) and 2L)/(j3, 4L), being 
in the numerator of the first and in the denominator of the second. As a result, 
these pairs of pairs are significantly antzcorrelated — one would expect a correlation 
coefficient of ~ —1/2, if all three raw-data points have roughly the same relative error 
- and they will thus tend to deviate from each other by more than would have been 
predicted from independent fluctuations with the given error bars. Furthermore, when 
L)/L is close to x*, 2L)/2Z/ is in turn not much smaller than £(/3,L)/L, so 
this anticorrelation acts on pairs of points having relatively nearby values of x. This 
could explain why the large contributions to x 2 come almost exclusively from the 
region x > 0.59, and why they are apparently completely random. Unfortunately, it 
seemed unfeasible for us to invert the large matrices (of order pa 70) that would be 
needed to take proper account of these correlations. Suffice it to say that if this is 
the correct explanation, then the observed large x 2 is simply spurious, and the fit is 
in reality good after all! Moreover, in this case the error bars on the extrapolated 
values £oo and Xst agg ,oo will be correct, as the method of "fake data sets" does take 
proper account of the aforementioned correlations. 

Modulo these caveats, therefore, we take as our preferred fit the one with n = 6, 
imin = 10 and L min = 128: we get 

1 + 2.234450e~ 1/:c + 80.120833e~ 2/x - 966.050470e" 3/x 

+ 10728.802555e- 4/:c - 49077. 175871e" 5/a ' + 73776.084137e~ 6/:r , 
(4-2) 

21 The L = 32 point at x 0.633 lies outside the range of x covered by the fit {x max « 0.630 
when L m i n = 128), so the negative deviation exhibited by this point may not be meaningful. 
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which is plotted in Figure |j. This fit is reliable in the interval x < x max ~ 0.629781. 
Using the value x* ~ 0.633888 derived from F^, we have 

F Xsta Jx*) = 3.172742 = 2 166573 (4.3) 

if we extrapolate ( |4.2| ) linearly for x > x max , or 

F Xstagg (x*) = 3.172345 = 2 L66555 (4.4) 

if we take (f4.2|) seriously also for x > x max . Either way, this is in excellent agreement 
with the prediction fSJ 7/1/ = 5/3 [cf. ([fsOp I- 



We can now compute the extrapolated values Xstagg, 00, using the functions F$ and 
F Xstagg given by our preferred fit as well as by several alternative fits that use more 
or less stringent choices of L min . (In all cases we take n = 5 for F^ and n = 6 for 
F Xstagg ] the results for larger n differ in almost all cases by < 0.6 standard deviations.) 
The statistical error bars on Xsta gg ,oo are computed as before by an auxiliary Monte 
Carlo process, and include errors of types (i) + (ii) + (iii). For each (3 we compute the 
weighted average ( |2.23[ ) of the various estimates Xstagg,oo, along with the statistical 



error bar ( |2.24j ) and the 1Z value ( |2.25|) . These are reported in Table |j. The statistical 
errors on Xstagg,™ are of order 2% (resp. 3%, 5%, 8%) at £«, « 1000 (resp. 2000, 5000, 
10000), that is, about twice as big as those on £oo. These extrapolated values from 
different lattice sizes at the same (3 are found to be consistent within statistical errors: 
only two of the 42 (3 values has an 1Z value too large at the 5% level; and summing 
all p values we have £7£ = 50.41 (75 DF, level = 98.7%). We don't know why this 
chi-squared is so small; it may be due in part to our overestimation of the raw-data 
error bars on £l, which leads to an overestimation of the errors of type (ii) on Xstagg,<x>- 
The discrepancies between the extrapolations with L min = 64, 128, 256 are again 
less than half the quoted statistical error in nearly all cases.0 We are thus reasonably 
confident that we have obtained quantitative control over the systematic errors due 
to corrections to scaling, and that their effect can be at most to double the quoted 
statistical errors. 



4.2 Analysis of Extrapolated Data 

In this section we analyze the behavior as (3 — > 00 of the extrapolated data for 
£ and Xstagg- in all cases we use the preferred extrapolations, namely the ones with 
i min = 10, L min = 128, and n = 5 (resp. 6) for £ (resp. Xstagg)- 

4.2.1 Correlation Length 

Our data are in clear agreement with the prediction of a critical point at (3 = 00. 
The correlation length ^ rises roughly like e 2/3 , and we initially thought that this 
was the exact asymptotic behavior. However, at /3 > 3.4 (£oo > 75), ^ begins to 

22 The extrapolations with L m i n — 32 frequently deviate from our preferred L m in = 128 extrap- 
olation by as much as w 1.2o\ These deviations are, once again, probably a correction-to-scaling 
effect on the borderline of statistical significance. 
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rise faster than this (Figure [|a), and this rise shows no sign of abating at least up to 
(3 ~ 5.9 (£oo ~ 15000). We therefore guessed a multiplicative logarithmic correction, 
i.e. £oo ~ e 213 (3 P for some power p > 0: see Figures ^p,c for p = 1/2 and p — 1, 
respectively. 

In order to distinguish between these scenarios, we need to make some assumption 
on the form of the additive corrections to the leading asymptotic behavior. Unfor- 
tunately we do not know how to carry out a low-temperature expansion around the 
(critical) zero-temperature state; but the simplest hypothesis is that there exists an 
expansion in powers of e _/3 , which corresponds to a minimum energy cost of one unit 
for an "overturned" spin. That is, we expect 

£«,(/?) = Ae 2f3 (3 p [l + aie"" + a 2 e" 2/3 + . . .] . (4.5) 

If we accept this Ansatz, a value p ~ 1 is clearly favored (Figure ||). A fit to the first 
two terms of ( (4.5[ ) with p — 1, using the data points with (3 > 2.95 (e~P < 0.052), 
yields A = 0.01814 ± 0.00006 and Aa x = 0.20051 ± 0.00225 (hence a x « 15) with 
X 2 = 13.10 (35 DF, level = 99.97%). 

On the other hand, Chris Henley (private communication) has suggested to us 
that the corrections to scaling might contain fractional powers of e~^\ 

£oo(£) = Ae 2/3 (3 P [l + a ie - Al/3 + a 2 e~ x ^ + ..] (4.6) 

with Xi < 1. Since e~ Xl/S ~ £ -Al//2 , such behavior would ordinarily arise from a 
correction-to-scaling exponent A = Ai/2; and our data for the finite-size-scaling func- 
tions show no evidence of a correction-to-scaling exponent anywhere near this small 
(see Section [4. 1 . 1| ) . Let us nevertheless consider the Ansatz ( f4.6| ) open-mindedly and 
see whether it can accommodate p = 0. In Figure ^a,b we plot CooW)/ e2lS versus e~ A/3 
for A = 1 and 0.5, respectively. With A = 1, the plot shows both strong curvature 
and a rather high slope near the origin; for (3 > 4.50 (e _/3 < 0.011) the data can be 
fit well by a straight line with A w 0.110 and a± ~ —15. With A = 0.5, the curvature 
and slope are less radical; the data for (3 > 4.50 (e~ 5/3 < 0.105) can be fit well by a 
straight line with A w 0.122 and a\ ~ —2.4. However, even this latter plot is nowhere 
near as convincing as Figure []. 

Finally, an anonymous referee has suggested to us that the corrections to scaling 
might be in powers of 1/(3: 



= Ae 2 ?{3 p 



(4.7) 



(see Section [74] for theoretical discussion). In Figure [8]a,b,c we plot £oo(/3)/(e 2/3 /3 p ) 
versus 1/(3 for p = 0, 1/2, 1. Clearly, our data — which lie in the range 0.15 ^ 1/(3 ^ 
0.4 — are so far from asymptotic that no conclusion can be drawn. In particular, the 
plot for p = (resp. p = 1) shows such a large negative (resp. positive) slope at the 
smallest available value of 1/(3 that the extrapolated intercept at 1/(3 = differs by 
nearly a factor of 2 from the last data point; and the plot for p = 1/2 shows a large 
third derivative (i.e. abrupt change from parabolic to flat) at 1/(3 ~ 0.2. 
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Let us note in passing that our data fit much less well the Ansatz ^ ~ exp(c/3 K ) 



used by some previous workers pq, 1C, 11]. Indeed, a log- log plot of log^ versus (3 



(Figure ||) shows significant curvature: the apparent exponent k, varies from « 1.69 
at small (3 to ~ 1.38 at larger (3. (This decrease is consistent with our conjecture 
that the true asymptotic value of the exponent k is 1.) Moreover, as noted in the 
Introduction, we have been unable to imagine any theoretical mechanism leading to 
£ ~ exp(c/5 K ) with k ^ 1. 

4.2.2 Staggered Susceptibility 



The staggered susceptibility is consistent with the believed exact behavior 
Xstagg,oo ~ £^ 3 , unmodified by any further powers of (3. To test this behavior quanti- 
tatively, we need to set error bars on the ratio Xstagg,oo/^ 3 - Since unfortunately we 
do not know the covariance between our estimates of £oo and Xsta gg ,oo, the best we 
can do is to use the triangle inequality to set an upper bound on the error bar for 
the ratio; this upper bound is of course a gross overestimate of the true error, since 
the estimates of £oo and Xsta gg ,<x are presumably strongly positively correlated. As 
a result, the error bars in all fits will be grossly overestimated, and the x 2 will be 
grossly underestimated; only the relative values of x 2 /DF have any significance. 

In Figure [U] we plot Xstagg,oo/CU 3 versus (3 (note the very narrow vertical scale); 
the error bars are those given by the triangle inequality, reduced by a factor of 10 
for visual clarity. We see that Xsta gg ,oo/ '£^ 3 varies slightly with f3, but appears to be 
tending to a constant ~ 2.67 as (3 — > oo. 

If we fit Xstagg,oc/ to the Ansatz A(3 P , the estimates of the power p are extremely 
small, of order 0.02, and statistically consistent with zero. This confirms our belief 
that there are no additional powers of (3 in the ratio Xsta gg ,oo/^ 3 - If we impose p = 
and fit Xstagg )00 /^ 3 to a constant A, we obtain A = 2.657 ± 0.051 (x 2 = 0.007, 13 
DF, level > 99.9%) using the data from (3 > 4.70. Of course, the error bar on A here 
is a gross overestimate, and the x 2 value is a gross underestimate. 

We can, of course, investigate directly the behavior of Xstagg,oo as a function of 
(3, without reference to £. This approach has the advantage that the error bars on 
Xstagg,oo are reliable. In Figure [11] we plot Xstagg,oo/(e ( - w/3 ^f3 q ) versus (3 for q = 0, §, §. 
The behavior is qualitatively similar to that observed in Figure [| for £, although the 
variation is somewhat sharper. If we try again the Ansatz 

Xstagg,oo(P) = Be^^" [l + b^ + b^ + . . .] , (4.8) 

a value q ~ 5/3 is favored (Figure |1^). The curvature is greater than in the cor- 
responding plot for £, but the linearity is still reasonable for e"' 3 < 0.03. A fit to 
( gg) with q = 5/3, using the data points with (3 > 3.60 (e^ < 0.0273), yields 
B = 0.00329 ±0.00003 and Bb x = 0.06661 ± 0.00124 (hence h w 20) with X 2 = 17.55 
(27 DF, level = 92%). We omit the plots based on a fractional-power additive correc- 
tion to scaling, which are similar to those shown for the correlation length (Figure 0). 
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4.3 Behavior of the Energy 

On theoretical grounds we expect that the infinite-volume energy per site has a 
low-temperature expansion of the form 

EM = Cl e-P + c 2 e- 2f5 + . . . . (4.9) 

Unfortunately we do not have access to E^, as we have made no attempt to ex- 
trapolate the energies (which are s/iort-distance quantities) to infinite volume. But 
examination of the finite-volume energies El indicates that for L > 512 (resp. 1024) 
the remaining L-dependence is less than w 2 x 10~ 4 (resp. 5 x 10~ 5 ); so the ener- 
gies appear to have almost reached their infinite-volume limits. In Figure [TBI we plot 
E/e" 13 versus e _/3 , using different symbols to represent different lattice sizes. A fit to 
the L = 1024 points with (3 > 5.00 {e~ p < 0.0067) yields c x = 0.21777 ± 0.00003 and 
c 2 = 1.65303 ± 0.00664 { X 2 = 5.09, 8 DF, level = 75%). The L = 1536 points are 
also compatible with this fit, and the L = 512 points differ only slightly from it. This 
provides good support for the Ansatz ( |4.9|) . 

We would like to make a warning concerning the use of reweighting methods^ 
|53|, 54, 55, 56] in which Monte Carlo runs at one temperature (3 are employed to 



generate data at another temperature j3' by reweighting with the factor exp[—(A(3)H], 
where A/3 = f3' — (3. This reweighting is of course always valid in principle; but one 
must be aware that the statistical error bars on the reweighted data grow rapidly as 
| A/3 1 grows, and the maximum \Af3\ for which one can obtain a not-too-large error 
bar gets smaller for larger L: 

!L~ d / 2 near a non-phase-transition point 
L~ l l u near a critical point (4.10) 
L~ d near a first-order phase-transition point 

Our Monte Carlo data illustrate this point in a very striking way. In all of our runs 
with L > 128 — totalling more than 60 million measurements — we did not observe 
even a single configuration (after the discard interval) with energy E = 0. In other 
words, none of our histograms for L > 128 — even those at our largest (3 value, namely 
(3 = 6.0 — have any overlap with the zero-temperature probability distribution. It 
follows that reweighting to zero temperature in these cases is nonsense. For L = 64, 
our runs at /3 — 4.50 and 5.00 (but not smaller (3) do show some configurations with 
E = 0, but the sample size of such configurations is very small: about 0.2% at /3 — 5.0. 
Reweighting to zero temperature in these cases would thus produce enormous error 
bars (at least if the error bars are computed correctly!). Only for L < 32 do we 
have a significant number of zero-energy configurations: for example, for L = 32 and 
(3 = 5.0, we found that 21.6% of the configurations have E — 0. Furthermore, even 
reweighting to nonzero temperatures is fraught with severe dangers. For example, 
although our run lengths are anywhere from 2 x 10 5 to 10 6 measurements, the energy 

23 These methods are sometimes called "histogram" (or "multiple- histogram" ) methods, but in 
fact the reweighting is most conveniently carried out without forming histograms! 
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histograms of a pair of runs typically show no overlap if |A/3| <; 0.3 (resp. 0.4, 0.5- 
0.6, 0.5-0.9, 1.2-1.4) for L = 1536 (resp. 1024, 512, 256, 128).Q It follows that 
reweighting beyond these limits is nonsense, and reweighting near these limits leads 
to huge statistical errors. 

5 Data Analysis: q = 3, Dynamic Quantities 

In this section we analyze the dynamic critical behavior of the WSK algorithm 
for the 3-state antiferromagnetic Potts model on the square lattice. 

5.1 Integrated Autocorrelation Times 

Examination of Table [I] indicates that the autocorrelation times T int ^2 and 
T%nt,e are bounded uniformly in (3 and L. (Indeed, their values are very small: 
T intM 2 < 5 and Tints < 4.) We conclude that critical slowing- down is completely 

' stagg 

eliminated. 

We can study the dynamic critical behavior in more detail by applying the stan- 
dard dynamic finite-size-scaling Ansatz 

T int>A ((3,L) « £(/3,L)^^(£(/?,L)/L) (5.1) 

to the observables A = M. 2 stagg and £. Here z intt A is a dynamic critical exponent, 
gA is an unknown scaling function, and (7a (0) = lim^ jo gA (x) is supposed to be finite 
and nonzero.0 Usually we would determine z intt A by plotting T intt A/£,{.L) Zint < A versus 
£(L)/L and adjusting z inty A until the points fall as closely as possible onto a single 
curve (with priority to the larger L values). But in our case the situation is much 
simpler: the dynamic critical exponents z intjM 2 and z intj£ are zero. 

In Figure |T3| we show the dynamic finite-size-scaling plot for T intjM 2 , assuming 
z int m 2 — 0. The data collapse is amazingly good, especially for a plot with no 

' stagg 

free parameters. Indeed, in all our Monte Carlo work on dynamic critical phenomena 
we have never observed a data collapse this good, even when we had the freedom to 
adjust z intyA . 

For Ti n t,e, by contrast, the data collapse is not so good: see Figure [L5|, where we 
again assume Zi ntt s = 0. Clearly there are huge corrections to dynamic finite-size- 
scaling for this observable. Even so, the points do appear to be converging as L — > oo 
to a limiting curve (which can be roughly traced using the L = 512 and L = 1024 
points) . 

24 The allowable | A/3| depends slightly on f3, getting smaller at smaller /3. For example, at L = 512 
the (3 = 5.7 run shows no overlap with the /3 = 5.1 run, which in turn shows no overlap with the 
P = 4.65 run, which in turn shows no overlap with the (3 = 4.2 run. 

25 We emphasize that the dynamic critical exponent Zi n t,A is in general different from the exponent 
z exp associated with the exponential autocorrelation time T exp |43L p7], pq] . 
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5.2 Autocorrelation Functions 



Now we want to test the more detailed dynamic finite-size-scaling Ansatz 

PA A(t;(3,L) « \t\- pA h A (t/r expA (P,L); £(0,L)/L) , (5.2) 

where pa is an unknown exponent and Ua is an unknown scaling function. If pa = 0, 
then (|5.2p can equivalently be written as 

p AA (t; (3, L) « h A (t/r mtA ((3, L) ; L)/L) . (5.3) 

In this latter situationP], t^a and T exPj A have the same dynamic critical exponent 
Zint,A = Zexp, and we furthermore have 



Tint,A 
T~exp,A 

where 



Fa L)/L) (5.4) 



F A (x) = lim j log/i A (t;x) . (5.5) 

Let us now test the Ansatz (|5~3| ) for the observable A = M^agg- (We restrict 
attention to this observable, since we already know that for A = £ the dynamic finite- 
size-scaling behavior is poor.) In Figure [16] we plot p M 2^ M ^ t (t) versus £/ T int,M 2 t , 
using all the data points. The points fall roughly on a single curve before falling into 
the statistical noise (which we expect to be of order ( n / r mt,A4^ agg ) where n is the 
run length, hence of order ±0.005). However, even at small £/ T mt,A^ t there are clear 
deviations from a single curve, indicating that the scaling function h M 2 depends 

stagg 

in a nontrivial way on its second argument £(L)/L. Therefore, in Figure |l7|a f we 
show the same plot with the data subdivided into "slices" of £(L)/L; the slices are 
chosen empirically so that the data points within a slice fall reasonably well onto a 
single curve modulo statistical noise. On each plot we also draw, for reference, a line 
corresponding to a pure exponential decay r int M 2 = T exPjM 2 . The data support 
the Ansatz ( |5.3j ) reasonably well, with each range of £ (L) / L defining roughly a single 
curve (until that curve falls into the statistical noise). The curves for small £(L)/L are 
close to straight (i.e. close to a pure exponential), while the curves for larger £(L)/L 
are increasingly convex.^] (Note that the rescaled horizontal axis ensures that the 
total area under each curve is 1. Therefore, the more convex curves must be below the 
straight curve for small t/r intM 2 but above it for large t/r intM 2 .) This means 

' stagg ' stagg 

that the ratio r int ^ /T e xp,M 2 t * s c l°se to 1 for small £(L)/L, and less than 1 for 
larger £(L)/L. It is conceivable that T int M 2^ lT e xp,M 2 t tends to 1 as £(L)/L — > 0; 
if true, this would mean that M 2 stagg truly becomes the "slowest mode" in the limit 
L^oo, f/L->0. 

26 Contrary to much belief, Zi n t a need not equal z exp . Indeed, if pa > 0, we have Zi n t A = 



(1 — pa) Zexp < z exp . See [43 [38| for further discussion. 

27 It is amusing to note that a similar behavior was observed in our study of the multi-grid Monte 
Carlo algorithm for the two-dimensional 0(3) cr-model |59[ |, We wonder whether it is a general 
phenomenon. 
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5.3 Relative Variance- Time Product 



Finally, let us look at the scaling fl2.26|) of the relative variance-time product. For 
each observable O — £, Xstagg we proceed as follows: For each run L) we form the 
relative variance (AOoo/Oqo) 2 on the extrapolated infinite-volume value coming from 
that run; we then multiply by L 2 x (^iterations — ^discard), which is a normalized 
measure of the CPU time invested in that run (after the discard interval); the result 
is, by definition, RVTP(/3,L). We then divide RVTP(/5,L) by ^((3) 2 [since d = 2 
and z intj0 = 0] and plot it versus ^ oo ( / 9)/L.0 The results are reported in Figure [l8| a 
for £ and Figure |l8| b for Xstagg- The scaling is reasonably good, though not perfect; 
this is not surprising, since our computed RVTP includes errors of types (i) + (ii) 
+ (iii) while the scaling formula ( |2.26| ) refers only to errors of types (i) + (ii).0Q 
Indeed, the fact that we see even modestly good scaling indicates that the errors of 
type (iii) are not dominant. 

We see that for £ the optimal value of ^oo/L is ~ 1, but the minimum is very 
flat: any value in the range 0.5 ^ ^oo/L ~ 10 is almost equally good. It is only for 
£oo/L < 0.3 that the RVTP rises sharply, by a factor of 10 or more. In other words, 
the only region in which one should not run is the region in which one traditionally 
always ran, namely ^oo/L ^1/6. For Xstagg the story is similar, but the minimum is 
somewhat sharper: the optimum is at C,oo/L « 0.4, and the RVTP rises by about a 
factor of 3 (resp. 10) as £oo/£ increases to ~ 10 (resp. decreases to ~ 0.1). 

6 Data Analysis: q = 4 

For q = 4 the story is very brief: simulations on L = 32, 64 agree within statistical 
error and show that £ ^ 2 uniformly as (3 — > oo (Table |]). Clearly there is no critical 
point in the physical region. Physically, there is so much entropy that the correlations 
decay exponentially even at zero temperature. This can be proven rigorously to occur 
on the square lattice for q > 7 || , and our simulations confirm Baxter's |23| prediction 
that it occurs in fact for q > 3. 

The autocorrelation times r int M 2 and T int s of the WSK algorithm are bounded 

' stagg 

uniformly in (3 and L, and indeed are almost constant: r intM 2 w 2.6 and Tints ~ 
3.5. 



28 In this latter computation, £oo(/3) is taken to be our best estimate (reported in Table [|), based 
on averaging all the runs at the given (3. 

29 The much cleaner graph shown in |2^, Figure 2] was computed by using the theoretical formulae 
for the propagation of errors under extrapolation |2l], ^2| , and includes only errors of types (i) + 
(ii). 

30 The anomalous L = 1536 point at £,00/ L w 16.6 corresponds to f3 = 6.0, where our statistics 
are poor and the extrapolation has almost broken down. Its relative error (33% for £ and 71% for 
Xstagg) is so large that both the extrapolated value and its error estimate are unreliable. 
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7 Discussion and Conclusions 



In this section we summarize our conclusions and discuss their theoretical impli- 
cations. We conclude by mentioning a few possible directions for future work. 

7.1 Behavior of the Correlation Length 

The numerical data presented in this paper (Section |4.2.1|) show clearly that the 
correlation length diverges as (3 —* oo approximately like £(/3) ~ e 2/3 . Since the 
fundamental variable in the Potts model is t = e J = e _/3 , an ordinary power- law 
critical point £ ~ (t—t c )~ u with t c — would correspond to £(/?) ~ e uf3 . Therefore, our 
result can be interpreted as indicating a power-law critical point at zero temperature 
with critical exponent v — 2. The fact that v is here a rational number reinforces 
our suspicion that this two-dimensional model can be solved exactly, at least in the 
sense of determining the exact asymptotic behavior as (3 — > oo. The exponent u = 2 
corresponds to an operator with scaling dimension X = 2 — l/V = 3/2, which is one 



of the possibilities proposed by Saleur |36], p. 248] — albeit not the one he considered 
most likely! 

On closer examination, however, the ratio £(/3)/e 2 ^ appears not to be asymptot- 
ically constant as (3 —* oo (see Figure ||a); rather, it begins to rise when /3 ~ 3.4 
(£ ~ 75), and this rise shows no sign of abating at least up to (3 ~ 5.9 (£ ~ 15000). 
Indeed, our data are compatible with an asymptotic behavior 

£(/3) = Ae 2/3 /? p [l + axe"' 3 + a 2 e" 2/3 + . . .] (7.1) 

with p « 1 (see Figure |B|). This corresponds to a power- law critical point with 
multiplicative logarithmic correction /5 P ~ | log(t — t c )\ p . The problem is to make 
theoretical sense of such a behavior. 

In the preliminary report of this work [13|], we asserted that a multiplicative log- 



arithmic correction (3 P ~ (\ogt) p with p integer (particularly p = 1) can occur in 
the renormalizat ion-group framework as a result of "resonance" between operators 
whose scaling dimensions are rationally related. This assertion is (we now realize) 
only half-true. For the susceptibility, specific heat and similar observables, it is in- 
deed true that multiplicative logarithmic corrections with positive integer powers p 
can occur as a result of resonance |B0|, ^T], However, we have been unable to 



devise any renormalization-group scenario in which the correlation length acquires 
such a multiplicative logarithmic correction in the absence of marginal operators. For 
example: 

1) Suppose that we hypothesize a scenario with one relevant variable t (eigenvalue 
A > 0) and one irrelevant variable u (eigenvalue —A), satisfying the RG flow equations 

^ = At + t 2 u (7.2a) 

du , , . 

- = —\u (7.2b) 
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where t(0) and u(0) are the couplings in the Hamiltonian, and t(l) and u(l) are the 
renormalized couplings after modes of momentum > e~ l have been integrated out. 
The solution is 



t{l) 
u{l) 



XI 



Ae 



-XI 



-Al 



(7.3a) 
(7.3b) 



where A = u(0). Let us now choose / so that t{l) = 1; this implies that e l = 
where £i is the correlation length at t — 1 and u = Ae~ xl ~ 0. Hence 



6 



l 



t(0) 



Alog(£/£i) 



l/A 



r^fi + o^iogO], 



(7.4) 



so that the presence of the irrelevant operator u (i.e. A^O) induces only an additive 
correction to scaling 0(t\ogt). 

2) Suppose, alternatively, that we hypothesize two relevant operators t, u with 
eigenvalues n\ and A, respectively, where n is an integer > 1 and A > 0, satisfying 
the RG flow equations 



dt 
Jl 
du 

~dl 



= nXt + u r - 



— = Xu 



The solution is 



t(l) = [t(0) +u(0) n l}e nXl 
u(l) = u{0)e xl 



(7.5a) 
(7.5b) 

(7.6a) 
(7.6b) 



Since there are two relevant operators, generically two couplings have to be adjusted 
in order to place the system at a critical point. But if some symmetry were to cause 
t(Q) to be exactly zero [or at least ^ u(0) n ], then a critical point can be reached 
by adjusting u(0) alone. So assume this, and let us analyze the resulting critical 
behavior. Let us choose I so that t(l) = 1; this implies that e l = where £i is the 
correlation length at t — 1 and u = l~ x l n ps 0. Hence 



nX 



so that 



— = u 1 I X | logu| 
?1 



-1/nA 



1 + O 



log log M 

log-u 



(7.7) 



(7.8) 



(One expects a further correction-to-scaling term 0( \ log it | 1 ' n ) arising from the fact 
that u(l) = l~ l l n 7^ 0.) Hence there is a multiplicative logarithmic correction, but its 
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power is negative; and there are very-slowly-decaying (logarithmic) additive correc- 
tions to scaling. 

3) In the presence of a marginally irrelevant operator u, multiplicative logarithmic 
corrections of either sign can be obtained.^ Suppose, for example, that the flow 
equations are 



dt 
Jl 



l(u)t 



du 

Ti = 

Then the solution is given by the implicit equation 

'«(0) du' 



I 



together with 



If we now assume that 



u(i) (3(u') 
u(0) 1 (u>) 



t(l) = t O fflp / ^rdu 
\Ju(l) P(U') 



(i) p(u>) 

A + 7!« + 7 2 m 2 + 7 3 u 3 + 
P2U 2 + P 3 u 3 + 



(7.9a) 
(7.9b) 

(7.10) 

(7.11) 

(7.12a) 
(7.12b) 



with A > and (3 2 > 0, and assume further that u(0) > is small enough so that 
(3{u) > for < u < u(0), we then find that 



t(l) 



u(l) 



const x t{0)e xl P l/l32 



1 + 



AsTi log I 
(31 I 



1 



Pi 1 



(7.13a) 
(7.13b) 



Setting now t(l) = 1 and e = as before, we obtain 



= r^liogtl-^/^ 



1 + 



/3s7i log 


log 






log* 





log* I 



(7.14) 



The exponent of the multiplicative logarithmic correction, —jx/^X, can thus be of 
either sign. Note, however, the presence of very-slowly-decaying additive corrections 
to scaling of the form O (log log tj log t) [with a universal coefficient] and 0(l/log£) 
[with a nonuniversal coefficient]. 

31 For a pedagogical discussion, see |63| , section 3.6]. Concrete manifestations of this phenomenon 
can be found in the four-dimensional N- vector model ]64| , |65[ , the three-dimensional tricritical 
TV-component model plfl , the three-dimensional iV-component ferromagnet with strong dipolar in- 
teractions ^3?J , and the two-dimensional 4-state Potts ferromagnet [[38| [39[ [70| . 
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Thus, we have been unable to devise any renormalizat ion-group scenario in which 
the correlation length acquires a multiplicative logarithmic correction with exponent 
p > in the absence of marginal operators; and we suspect that no such scenario 
exists. On the other hand, if marginal operators were present in our model (as in 
scenario #3), then one of their effects would be to induce 1/ logL corrections to scaling 
in the finite-size-scaling functions^, and we see no evidence of corrections decaying 
anywhere near so slowly. On the contrary, the corrections to scaling are here almost 
undetectable, and they appear to decay like L~ A with A > 1 (see Section [4.1. 1| ). 

Finally, an anonymous referee has pointed out to us the possibility that the im- 
portant microscopic variable is not t = e _/3 but rather t = 1/(3, and that this t is a 
marginally relevant operator: 

^ = B(t) = b 2 t 2 + b 3 t 3 + ht 4 + ... (7.15) 

with 62 > 0, i.e. a situation of asymptotic freedom. (This sounds implausible at first 
sight for a discrete-spin model, but is not impossible.) The solution to (|7.15|) is given 
by 

/■*(Q dt' 

I = / JL- . (7.16) 
Jt(p) B(t>) v ; 

If we set t(l) = 1 and t(0) = 1/(3 and use e l = as before, we obtain 



f = Ae^ (3- b3/bl 



1 + 0(1 



(7.17) 



where A is a nonperturbative constant, and the corrections in inverse powers of (3 
can be computed from the coefficients 64, 65, . . . . On the other hand, the corrections 
to finite-size scaling are determined (as usual) by irrelevant operators, and so decay 
as inverse powers of L (provided there are no marginally irrelevant operators). So 
this scenario, unlike scenario #3, is not ruled out by our finite-size-scaling data (Sec- 
tion [4.1. lp . However, if (|7.17[ ) is indeed the true behavior, then our data are so far 
from asymptotic that they give no useful information about A and p = —bj > /b\ (see 
Section [4.2. 1| and in particular Figure El). So we cannot rule out this scenario, but 



neither can we obtain any evidence in its favor. 

In summary, we really do not understand the theoretical basis for an asymptotic 
behavior of the form ( [7.1|) with p > 0. 

An alternative possibility is that the true asymptotic behavior is £(/3) ~ e 2/3 
without multiplicative logarithmic corrections. This could happen in either of two 
ways: 

a) The rise seen in Figure^a at (3 ^ 3.4 might be spurious, i.e. an artifact of some 
undetected systematic error in our extrapolation method. Indeed, it is suspicious 
that this rise begins at roughly the same correlation length (£ w 75) where our 

32 See 0, Section 3] for a detailed theoretical study of the 1 / log L corrections in the finite-size- 
scaling functions when a marginally irrelevant operator is present; and see |70| , Figure 6] for an 
illustration of their practical effect in the two-dimensional 4-state Potts ferromagnet. 
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extrapolation method begins to play a central role.0 To test whether this rise is real, 
we produced a traditional finite-size-scaling plot in which £ / e 2/3 is plotted versus £/L 
(Figure [19]). One sees clearly that the points do not fall on a single curve, and that it 
is the largest lattices that deviate the most from the others (see particularly the range 
0.25 ^ £/L ^ 0.4); the data are most definitely not compatible with convergence to 
a limiting FSS curve. Rather, the deviations reflect precisely the rise of £/e 2 ^ when 
£ i5 100. We therefore think that the rise observed in Figure |^a is real. 

b) The rise seen in Figure |3|a might level off at some ft > 6. This is perfectly 
possible, but it would mean that the corrections to the leading asymptotic behavior 
either have unusually strong amplitude or else decay more slowly than the hypothe- 
sized e _/3 . If the corrections behave as in ( |7.1|) , the coefficient ai would have to be 
approximately —15 in order to obtain a decent fit between ft = 4.5 (e~@ ~ 0.011) and 
ft = 6.0 (e _/3 ~ 0.002), and the coefficients a 2 and a 3 would have to be large as well 
(see Figure 0a). On the other hand, if we allow additive corrections to scaling that 
are fractional powers of e _/3 , 

£(J3) = Ae 2p [l + aie~ Al/3 + a 2 e" A2/3 + . . .] , (7.18) 

then the data for ft > 4.5 can be fit reasonably well with Ai ~ 0.5 and a\ « —2.4 
(see Figure |7|b). Now such nonanalytic corrections to scaling can arise routinely 
from irrelevant operators: in the case at hand, an additive correction e~ Xl/S ~ £~ Al / 2 
corresponds to a correction-to-scaling exponent A = Ai/2 « 1/4. Unfortunately, 
our study of the finite-size-scaling function (Section |4.1.1|) gives no indication of any 
correction to scaling with A < 1. So it is highly unlikely that a correction e~ Al/3 
with Ai ~ 0.5 could arise from this mechanism. We do not know whether some other 
mechanism might lead to such behavior. 

In conclusion, two distinct Ansatze on the large-/? asymptotic behavior of the 
correlation length — ( 7.1j) with p « 1, and ( |7.18j ) with Ai ~ 0.5 — are compatible 



with our data, but both present difficulties of theoretical interpretation. We hope 
that someone will be able to sort this out, and that the numerical results presented 
here will serve as useful clues toward the exact solution of this model, possibly with 
the help of the methods of conformal field theory |5|, [|. Our estimate of the universal 
quantity x* ~ 0.633888 (see also f54| ) provides another constraint in determining the 
universality class of this model. 



7.2 Prospects for Future Work 

Here are some possible directions for future work on the 3-state square-lattice 
Potts antiferromagnet: 

1) Study small lattices (e.g. L — 4, 8, 16, 32, 64) with very high statistics in order 
to obtain quantitative information on the corrections to scaling. Salas and Sokal j34] 
have recently done this at ft = oo (corresponding to £(L)/L = x* ~ 0.634), but 

33 From Figure |l| we see that the finite-size corrections are negligible when £/L < 1/6. Inspection 
of Table [l] shows that we have raw data satisfying £/L < 1/6 up to (3 — 3.50 (£ « 94), but not 
beyond that. 
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the corrections to scaling at this value of £(L)/L are quite weak. As can be seen 
in Figures 0(b) and |3|(b), the corrections are much stronger in the interval 0.54 < 
£(L)/L ^ 0.62, so it would be useful to get higher statistics in this region. 

2) Study large lattices (e.g. L = 1024, 1536,2048) at f/L < 1/6 (where the finite- 
size effects are negligible, see Figure [3]) in order to verify that the observed rise of 
£(/?)/e 2/3 at f3 > 3.4 (£ > 75) is real and not merely an artifact of our extrapolation 
method. Using L = 2048, we could hope to reach /3 ~ 4.1 (£ ~ 320) with "essentially 
infinite- volume" simulations. By this time the rise of £,{/3)/e 213 is about 4% (see 
Figure ||a) and so should be clearly detectable. 

Future studies of this model should also correct two defects in the present work: 

3) We measured here the susceptibility and correlation length associated to the 
staggered magnetization A4 s t a gg, which is the most relevant operator in this model 
(Vstagg = 1/3) and hence has the most strongly divergent susceptibility (jstagg/v = 
2 — Vstagg = 5/3). What we failed to notice (until the runs had already been made and 
it was therefore too late!) is that the uniform magnetization M. u is also a relevant 
operator (jj u = 4/3 |3(| |33|, [34|) with a divergent susceptibility (y u /v = 2/3). Future 
studies should measure it as well. Finally, the staggered polarization V s tagg (see |3B|, [34] 
for the definition) is also a relevant operator {f]-p atagg = 3) albeit with a non-divergent 
susceptibility {y-p sta g /v = —1). By measuring its correlation function at several 
different low momenta, it ought to be possible to check the prediction rjp gtagg = 3 
despite the non-divergence of the corresponding susceptibility. 

4) In the present paper we didn't bother to measure the cross-correlations between 
M stagg and T s tagg'i as a result, we were unable to assign statistically valid error bars to 
£ (instead we used the triangle inequality to obtain an overly conservative error bar). 
This was a serious mistake, as it prevented us from distinguishing clearly between 
statistical fluctuations and systematic errors (arising from corrections to scaling or 
other sources): see Section |4.1.1| . Measuring the cross- correlations would allow one 
to determine the correct error bars not only on £ but also on combinations such as 
x/£ 5 / 3 ; this could lead to a very sensitive test of the conjecture that x/£ 5//3 — > const 
as (3 -> oo (Section pLT2p . 

A Proof of a Correlation Inequality for Antiferro- 
magnetic Potts Models on a Bipartite Graph 

Let V be a finite set of sites; we shall consider a g-state Potts model on V consisting 
of spins <Tj G {1, . . . , q} for i 6 V, interacting via a Hamiltonian 

n = -E-^Wj • (a.i) 

m 

Here the sum runs over all pairs i, j G V (each pair counted once), and {Jij}i,jev is 
some specified set of couplings. 

In this appendix we shall prove some correlation inequalities for the Potts model 
( [A.I]) . Our technique will be to embed a field of Ising spins into the given Potts 
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model — using, in fact, precisely the Wang-Swendsen-Kotecky (WSK) embedding 
discussed in Section |2]3| — and to exploit the well-known Griffiths inequalities Jn], [72| 
for the induced Ising model. It is amusing that the WSK embedding can be used 
both as a Monte Carlo algorithm and as an analytical technique. Similar proofs of 
correlation inequalities based on the embedding of Ising or XY variables can be found 
0,0, 0,1761, ]77J. 



m 



We begin by introducing a redundant (but perfectly legitimate) parametrization 
of the g-state Potts spin Oi in terms of a g-state Potts spin oji and an Ising spin Tj. 
Let u)i be uniformly distributed in {1, . . . , q}, and let Tj be uniformly distributed in 
{— 1, +1}; we then define 

a, = (H* if -. = lor2 (A.2) 
[ uji if LUi > 3 

It is easy to see that <Tj is uniformly distributed in {1, . . . , q}: indeed, each of the 
q possible values of <Tj arises from exactly two of the 2q possible values of Tj). 
Performing this construction independently at each site i G V, we construct the a 
priori measure for the spins {ui, Ti} iG v, and thus also for the spins <7j defined by 
( |A.2j ) as functions of (u>i, Tj). The desired probability distribution is then obtained by 
multiplying this a priori measure by Z^ 1 exp[—H.({a})}, where 7i is given by ( |A.1| ) 
and {o} is defined by ( |A.2j ). 

Let us first compute the probability distribution of the set of spins {r} conditioned 
on the set of spins {o}. Here the Boltzmann weight factor Z~ l exp[— 7i({<j})] plays 
no role, since the {o} are being held fixed; the conditional distribution is the same 
as in the a priori measure. We thus have, independently for each site i, 



2(jj — 3 if <Tj = 1 or 2 

±1 with equal probability if Oi > 3 



(A.3) 



It follows that the conditional expectation of TjTj given {a} is 

E{nTj\{<r}) = I{a i <2)I{a j <2){28 aua .-l) (A.4a) 

+ 1 if o"j = Oj = 1 or Oi = Oj = 2 
-1 if Oi = 1, Oj = 2 or Oi = 2, Oj = 1 (A.4b) 
in all other cases 

where /(• • •) denotes the indicator function of the specified event. Next we want to 
take the unconditional expectation of (A. 4). Since the probability distribution of the 
{o} is invariant under global permutations of {1, . . . , q}, we have in particular that 
the joint probability distribution of o^ and Oj is given by 

Prob(^) = 1 5^ + (A.5) 



where p = (5 CT . iCr .). Averaging (A. 4) over the distribution ( |A.5|) , we find that 



v 3 ' g(g-i) q\ g-i K ' 
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Note that on the right-hand side of this identity we have precisely the two-point 
correlation function of our Potts model [cf. (|2.3|) 1, multiplied by 2/q. 

On the other hand let us compute the probability distribution of the set of spins 
{r} conditioned on the set of spins {ui}. Note first that 

= I(ui > 3) Ifa > 3) 8 Ui ^ + I(ui < 2) I(uj < 2) 1 + ^ iTj . (A.7) 

Inserting this into (A.l), we see that the model of spins {r} conditioned on {a;} is an 
Ising model with interactions 

Jtf = ij ij I(u i <2)I(u j <2). (A.8) 

(The factor ~ comes simply from the difference between the conventional Ising and 
Potts normalizations of couplings.) Below, we shall apply various correlation inequal- 
ities to this induced Ising model. 

Let us consider some special cases: 

Example 1 (ferromagnetic Potts model). Assume that > for all i,j G V. 
Then the couplings J*f are also > 0. So Griffiths' first inequality |TT| applied to the 
induced Ising model ( |A.8| ) implies that 

E(nT j \{u}) > for all {cu} . (A.9) 

In particular, averaging this over {u} we deduce that (T»Tj) > 0; and hence by the 
identity ( |A.6|) we have 

( ^ 5 ^Z~ > 0. (A.10) 

This is a "Griffiths' first inequality for ferromagnetic Potts models" . 

Of course, this inequality can be proven somewhat more simply using the Fortuin- 
Kasteleyn representation: the left-hand side of (|A.10|) is equal to a connection prob- 
ability in the Fortuin-Kasteleyn bond variables, and thus is manifestly nonnegative. 
Indeed, this proof is valid for all real q > 0. So nothing much is gained by WSK 
embedding in this case; the real value of the method arises in the next case: 

Example 2 (antiferromagnetic Potts model on a bipartite graph). Assume that 
the set of sites V can be partitioned as V = A U B in such a way that 

{> if i,jeA 
>0 iii,jeB (A.ll) 
< if i G A, j e B 

(In particular, the pure antiferromagnet would have Jjj = for i,j G A and for 
i,j G B.) Now define 



-Ti 
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It follows that the model of spins {r'} conditioned on {u} is an Ising model with 
interactions 



+1 HiJeA 
-1 iii,jeB 
-1 if i e A, j e B 



J-f= \JijI{ui < 2)I(uj < 2) x { +1 if i,j e B } > 0. (A.13) 



So Griffiths' first inequality |7TJ applied to this ferromagnetic Ising model implies that 

E{tIt^{uj}) > for all {u} . (A. 14) 

Averaging over {a;} we deduce that (r/rj) > 0, and hence by the identity (|A.6|) that 

, , N r>0 if i,je A 

A Utw) >0 if,./.: I'> (A.15) 



q - 1 2 



< if i e A, j e B 



This is a kind of "Griffiths' first inequality for antiferro magnetic Potts models". 
In particular, for antiferromagnetic Potts models on the square lattice, (|A.15 ) 



implies that the two-point correlation function G(x, y) has sign (—ly x ~ y \ where \x— y\ 
denotes the i 1 norm. One consequence of this is that the Fourier-transformed two- 
point function G(p) satisfies 

\G(p)\ < G((-k,tt)). (A.16) 

Example 3 (comparison-to-Ising inequality). Let us return to the general case 
of arbitrary { J^}. From (|A.8| ) it follows immediately that 

\J?f\ < l\Jij\ , (A-17) 
so by Griffiths' comparison inequality JF^| we have 

lE^TiTjliu})] < (eiej)ising,{\j/2\} for all {w} , (A.18) 

where ( • • • ) ising, {\J/2\} denotes the expectation in an Ising model with couplings 
{|| Jij\}. Averaging over {a;} and using ( |A.6| ), we conclude that 

^2 I j ^ | {^j)lsing,{\J/2\} ■ (A.19) 

Thus, the two-point function in an arbitrary Potts model can be bounded above by 
the corresponding two-point function in a ferromagnetic Ising model.0 

Of course, for a ferromagnetic Potts model a much stronger result is true, namely that each 
correlation function G(x,y;q) is a decreasing function of q (for q > 1) at fixed {Jy}. This fol- 
lows from the Fortuin-Kasteleyn representation combined with the FKG inequality for the random- 



cluster model j78|, |79|]. So the real interest of (A.19) is for the antiferromagnetic or mixed ferromag- 
netic/ antiferromagnetic Potts models. 
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For example, one consequence of ( A.19 ) is that the transition temperature of the 
q = 3 Potts antiferromagnet on the triangular lattice must satisfy 



Ptrans (TRI g = 3 Potts AF) > 2/3 crit (TRI Ising ferro) = -log 3 = 0.549306... . 

(A.20) 

This bound is of course satisfied by the numerical estimate /5 traris (TRI q = 3 Potts AF) 
1.594 [ 30 1 . Likewise, for the q = 3,4 Potts antiferromagnets on the simple-cubic lat- 
tice, one must have 



Ptrans(SCq = 3,4 Potts AF) > 2/3 cHt (SC Ising ferro) 



(A.21) 



which is satisfied by the numerical estimates /3t rcm s(SC q = 3 Potts AF) « 0.816 ]TT 
|82|, |8|, g, /3 tmns (SC g = 4 Potts AF) « 1.43 |85j, and /? crii (SC Ising ferro) r 



0.222 86 



B Proof of Ergodicity of the WSK Algorithm at 
T = on a Bipartite Graph 

Let G = (V, E) be a finite undirected graph with vertex set V and edge set E. 
Then G is said to be bipartite if the vertex set V can be partitioned as V = A U B in 
such a way that every edge e G has one endpoint in A and the other endpoint in 
B (i.e. there are no A— A or B-B edges). 

Example 1. Any finite subset of the simple (hyper-)cubic lattice Z d , with free 
boundary conditions, defines a bipartite graph. 

Example 2. A box of size L\ x L 2 x . . . x in Z d , with periodic boundary 
conditions, defines a bipartite graph if and only if all the side lengths L 1; L 2 , . . . , 
are even. 

The goal of this Appendix is to prove the following theorem: 

Theorem B.l Let G be a bipartite finite undirected graph, and let q be an integer > 
2. Then the Wang-Swendsen-Kotecky algorithm for the q-state Potts antiferromagnet 
on G at zero temperature (i.e. for q-colorings of G) is ergodic. 



PROOF. We will prove, by induction on q, that the WSK algorithm for g-colorings 
(hereafter called WSK-g) is ergodic on G and on all its subgraphs. 

Case q = 2. Since G is bipartite, so are all its subgraphs H. The WSK algo- 
rithm acts independently on each connected component of H, so it suffices to prove 
the ergodicity for each connected component. But a connected bipartite graph has 
precisely two 2-colorings, which are related by a global interchange of the two colors; 
and this global interchange can trivially be realized by a WSK move. 

Inductive step. Let q > 3, and suppose that WSK-(g — 1) is ergodic on G and all 
its subgraphs; we shall prove the same for WSK-g. So let H = (V, E 1 ) be a subgraph 



33 



of G, and define A' — An V, B' = B H V. Define the reference configuration to 
be the g-coloring of H in which all sites in A' are colored 1 and all sites in B' are 
colored 2. Let {o~ x } xe yi be an arbitrary g-coloring of H, which we call the target 
configuration. We shall show that the reference configuration can be transformed 
into the target configuration by a finite sequence of WSK-g moves. (This is sufficient 
to prove ergodicity, since the inverse of a WSK move is also a WSK move.) 

(a) Step 1. Choose l,q as the pair of colors to be used in the WSK move, and 
focus attention on all sites x G A' such that the target configuration has a x = q. 
All these sites are of course currently colored 1, and all their neighbors (which are 
in B') are currently colored 2. So each of these sites is a singleton (i.e. a one-site 
connected component) in the subgraph of H formed by those sites currently colored 1 
or q. Therefore, each of these sites can be recolored q by a WSK move, while leaving 
all other sites unchanged. 

(b) Step 2. Choose 2, q as the pair of colors to be used in the WSK move, and 
focus attention on all sites x G B' such that the target configuration has a x = q. All 
these sites are of course currently colored 2; and all their neighbors (which are in A' 
and have target color ^ q) are currently colored 1. So each of these sites is a singleton 
in the subgraph of H formed by those sites currently colored 2 or q. Therefore, each of 
these sites can be recolored q by a WSK move, while leaving all other sites unchanged. 

(c) Step 3. All sites with target color q are now colored q, while the remaining 
sites are now colored 1 or 2. The latter sites (along with their corresponding edges) 
define a subgraph K C H; and by the inductive hypothesis they can be given their 
target colors (which lie in {1,2, ... ,q — 1}) by a sequence of WSK-(g — 1) moves 
(which are of course also WSK-g moves). I 

Note. After completion of this work, we learned that the same theorem (as well 
as some generalizations of it) was obtained independently by Burton and Henley [[33 
Their proof is essentially the same as ours. 
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Table 1: Our Monte Carlo data for the 3-state Potts antiferromagnet. "Total" is 
the total number of WSK iterations performed; "Discard" is the number of itera- 
tions discarded for equilibration. Error bars (one standard deviation) are shown in 
parentheses. [First page of 3-page table] 



T 

Li 


P 


lotai 


Discard 




Xstagg 


E 


T int,M 2 Btagg 




128 


4.10 


200000 


10000 


75 


419 


(0.204) 


2594.69 


( 3 


94) 


0.00406992 


(0.00000280) 


4.32 


(0.10) 


2 


69 


(0.05) 


128 


4.20 


200000 


10000 


76 


508 


(0.205) 


2636.64 


( 3 


92) 


0.00363639 


(0.00000262) 


4.52 


(0.11) 


2 


64 


(0.05) 


128 


4.30 


1000000 


10000 


77 


427 


(0.090) 


2670.89 


[ 1 


64) 


0.00325960 


(0.00000107) 


4.29 


(0.04) 


2 


62 


(0.02) 


128 


4.40 


1000000 


10000 


78 


050 


(0.089) 


2697.29 


( 1 


61) 


0.00291971 


(0.00000101) 


4.27 


(0.04) 


2 


61 


(0.02) 


128 


4.50 


1000000 


10000 


78 


772 


(0.090) 


2723.43 


( 1 


59) 


0.00262197 


(0.00000095) 


4.32 


(0.04) 


2 


59 


(0.02) 


128 


4.60 


1000000 


10000 


79 


136 


(0.090) 


2742.07 


( 1 


59) 


0.00235324 


(0.00000090) 


4.37 


(0.05) 


2 


59 


(0.02) 


128 


4.70 


1000000 


10000 


79 


505 


(0.091) 


2756.45 


( 1 


58) 


0.00211402 


(0.00000084) 


4.41 


(0.05) 


2 


56 


(0.02) 


128 


4.80 


1000000 


10000 


79 


952 


(0.090) 


2771.60 


( 1 


56) 


0.00190149 


(0.00000080) 


4.40 


(0.05) 


2 


56 


(0.02) 


128 


4.90 


1000000 


10000 


80 


091 


(0.091) 


2781.79 


( 1 


56) 


0.00171139 


(0.00000075) 


4.42 


(0.05) 


2 


56 


(0.02) 


128 


5.00 


1000000 


10000 


80 


174 


(0.091) 


2788.86 


( 1 


56) 


0.00154123 


(0.00000071) 


4.42 


(0.05) 


2 


53 


(0.02) 


128 


5.10 


1000000 


10000 


80 


443 


(0.090) 


2797.74 


( 1 


54) 


0.00138704 


(0.00000067) 


4.37 


(0.05) 


2 


52 


(0.02) 


128 


5.20 


1000000 


10000 


80 


612 


(0.091) 


2805.76 


( 1 


54) 


0.00125012 


(0.00000064) 


4.40 


(0.05) 


2 


55 


(0.02) 


256 


2.50 


200000 


10000 


13 


104 


(0.422) 


201.34 


( o 


81) 


0.03295828 


(0.00000517) 


1.57 


(0.02) 


3 


34 


(0.07) 


256 


2.60 


200000 


10000 


16 


185 


(0.356) 


273.13 


[ 1 


08) 


0.02839657 


(0.00000461) 


1.56 


(0.02) 


3 


18 


(0.07) 


256 


2.70 


200000 


10000 


19 


126 


(0.325) 


369.05 


( 1 


51) 


0.02450313 


(0.00000423) 


1.69 


(0.03) 


3 


15 


(0.06) 


256 


2.80 


200000 


10000 


23 


045 


(0.293) 


502.20 


( 2 


08) 


0.02119791 


(0.00000393) 


1.76 


(0.03) 


3 


16 


(0.06) 


256 


2.90 


200000 


10000 


28 


119 


(0.269) 


693.32 


( 2 


93) 


0.01837763 


(0.00000362) 


1.90 


(0.03) 


3 


16 


(0.06) 


256 


3.00 


200000 


10000 


34 


305 


(0.258) 


954.62 


( 4 


09) 


0.01599052 


(0.00000333) 


2.09 


(0.04) 


3 


13 


(0.06) 


256 


3.10 


200000 


10000 


41 


701 


(0.260) 


1313.21 


( 5 


88) 


0.01394816 


(0.00000305) 


2.48 


(0.04) 


3 


10 


(0.06) 


256 


3.20 


200000 


10000 


50 


453 


(0.269) 


1803.02 


( 8 


21) 


0.01219458 


(0.00000281) 


2.96 


(0.06) 


3 


05 


(0.06) 


256 


3.30 


200000 


10000 


60 


839 


(0.286) 


2439.41 


10 


78) 


0.01069957 


(0.00000263) 


3.35 


(0.07) 


3 


08 


(0.06) 


256 


3.40 


200000 


10000 


72 


254 


(0.316) 


3198.35 


13 


71) 


0.00940398 


(0.00000238) 


3.99 


(0.09) 


2 


94 


(0.06) 


256 


3.50 


1000000 


10000 


84 


685 


(0.150) 


4036.07 


' 6 


91) 


0.00829150 


(0.00000096) 


4.38 


(0.05) 


2 


92 


(0.03) 


256 


3.60 


200000 


10000 


96 


720 


(0.365) 


4854.04 


16 


67) 


0.00732740 


(0.00000203) 


4.48 


(0.11) 


2 


85 


(0.06) 


256 


3.70 


200000 


10000 


107 


817 


(0.380) 


5595.86 


16 


87) 


0.00648727 


(0.00000188) 


4.59 


(0.11) 


2 


82 


(0.05) 


256 


3.80 


200000 


10000 


117 


698 


(0.392) 


6228.98 


16 


33) 


0.00576244 


(0.00000175) 


4.50 


(0.11) 


2 


81 


(0.05) 


256 


3.90 


200000 


10000 


125 


908 


(0.414) 


6756.45 


16 


21) 


0.00512698 


(0.00000163) 


4.77 


(0.12) 


2 


78 


(0.05) 


256 


4.00 


200000 


10000 


132 


919 


(0.413) 


7173.99 


15 


25) 


0.00456912 


(0.00000151) 


4.62 


(0.11) 


2 


72 


(0.05) 


256 


4.10 


200000 


10000 


138 


468 


(0.409) 


7524.03 


14 


19) 


0.00407866 


(0.00000140) 


4.40 


(0.11) 


2 


64 


(0.05) 


256 


4.20 


200000 


10000 


143 


448 


(0.409) 


7812.70 


13 


50) 


0.00364854 


(0.00000133) 


4.30 


(0.10) 


2 


70 


(0.05) 


256 


4.30 


200000 


10000 


146 


776 


(0.404) 


8024.93 


12 


99) 


0.00326371 


(0.00000124) 


4.29 


(0.10) 


2 


66 


(0.05) 


256 


4.40 


200000 


10000 


150 


305 


(0.412) 


8207.19 


12 


71) 


0.00292610 


(0.00000116) 


4.36 


(0.10) 


2 


58 


(0.05) 


256 


4.50 


1000000 


10000 


152 


769 


(0.180) 


8361.37 


( 5 


39) 


0.00262364 


(0.00000048) 


4.33 


(0.04) 


2 


62 


(0.02) 


256 


4.60 


1000000 


10000 


154 


560 


(0.181) 


8472.30 


' 5 


32) 


0.00235575 


(0.00000045) 


4.39 


(0.05) 


2 


56 


(0.02) 


256 


4.65 


1000000 


10000 


154 


914 


(0.181) 


8512.61 


( 5 


26) 


0.00223281 


(0.00000043) 


4.36 


(0.05) 


2 


57 


(0.02) 


256 


4.70 


1000000 


10000 


155 


773 


(0.179) 


8560.34 


( 5 


21) 


0.00211659 


(0.00000042) 


4.35 


(0.05) 


2 


57 


(0.02) 


256 


4.80 


1000000 


10000 


157 


217 


(0.179) 


8645.33 


( 5 


05) 


0.00190282 


(0.00000040) 


4.25 


(0.04) 


2 


54 


(0.02) 


256 


4.90 


1000000 


10000 


157 


807 


(0.179) 


8696.27 


( 5 


02) 


0.00171308 


(0.00000038) 


4.27 


(0.04) 


2 


53 


(0.02) 


256 


5.00 


1000000 


10000 


158 


676 


(0.180) 


8750.44 


( 5 


01) 


0.00154114 


(0.00000036) 


4.35 


(0.05) 


2 


54 


(0.02) 


256 


5.10 


1000000 


10000 


159 


453 


(0.182) 


8796.18 


( 5 


02) 


0.00138804 


(0.00000034) 


4.43 


(0.05) 


2 


57 


(0.02) 


256 


5.20 


1000000 


10000 


160 


061 


(0.182) 


8833.00 


( 4 


98) 


0.00125104 


(0.00000032) 


4.42 


(0.05) 


2 


56 


(0.02) 


256 


5.30 


1000000 


10000 


160 


395 


(0.182) 


8864.52 


( 4 


96) 


0.00112784 


(0.00000030) 


4.43 


(0.05) 


2 


53 


(0.02) 


256 


5.40 


1000000 


10000 


160 


835 


(0.182) 


8892.90 


( 4 


93) 


0.00101687 


(0.00000029) 


4.40 


(0.05) 


2 


52 


(0.02) 


512 


2.80 


500000 


10000 


23.428 


(0.577) 


505.57 


( 1-25) 


0.02119445 


(0.00000123) 


1.54 


(0.01) 


3.22 


(0.04) 


512 


2.90 


500000 


10000 


27.679 


(0.508) 


690.13 


( 1-74) 


0.01838339 


(0.00000113) 


1.61 


(0.01) 


3.17 


(0.04) 


512 


3.00 


500000 


10000 


33.500 


(0.442) 


949.08 


( 2 - 41 ) 


0.01598973 


(0.00000104) 


1.64 


(0.02) 


3.14 


(0.04) 


512 


3.10 


200000 


10000 


41.719 


(0.619) 


1325.44 


( 5.48) 


0.01394838 


(0.00000152) 


1.73 


(0.03) 


3.07 


(0.06) 


512 


3.20 


200000 


10000 


50.521 


(0.559) 


1830.24 


( 7.59) 


0.01220209 


(0.00000139) 


1.81 


(0.03) 


2.99 


(0.06) 


512 


3.30 


500000 


10000 


61.816 


(0.330) 


2552.67 


( 6.83) 


0.01070391 


(0.00000080) 


2.03 


(0.02) 


3.00 


(0.04) 


512 


3.40 


500000 


10000 


75.572 


(0.320) 


3563.37 


( 9.72) 


0.00941670 


(0.00000075) 


2.29 


(0.02) 


2.99 


(0.04) 


512 


3.50 


500000 


10000 


92.478 


(0.330) 


4976.76 


(14.11) 


0.00830439 


(0.00000068) 


2.76 


(0.03) 


2.90 


(0.04) 



Table 1: 



[Second page 

40 



of 3-page table] 



L 


P 


Total 


Discard 




Xstagg 




T int,M 2 st 




7~int,£ 


512 


o 
6 


60 


500000 


10000 


1 1 o 

112 


O A C 

o45 


/n 
v 


O C A\ 

o50) 


0836 


C O 

58 


/ in 

iy 


07\ 

o7J 





AATO A A d A 

007o4161 


( A 

(U 


A AAA A C\CZ O \ 
0000006O ) 




3 


O^ 

26 


/ A 

v 


04) 



2 


OO 

82 


I C\ AO\ 

(O.Oo) 


512 


o 



a tr 
DO 


500000 


10000 


1 oo 

12o 


ooo 
2oo 


( n 
v 


ODD j 


TOO O 

7933 


no 
08 


( oo 
22 


46) 


A 




oo6yoy tb 


Cn 
(0 


A A A A A A/? 1 ^ 

00000061 J 


'_> 

3 


A 

60 


( A 

v 


05) 



2 


rr 
85 


(O.Oo) 


512 


O 
O 




500000 


10000 


loo 


iyi 


( n 

V) 


o82 ) 


yioo 


68 


( OK 

2o 


Q A \ 

o4 ) 


A 




nnft Knxnft 
UU65U5U6 


/"n 
v u 


0UUU0U58 ) 


6 


8 / 


( (\ 

v 


05) 



z 


1 
81 


/ A C\0\ 

(U.Oo ) 


512 


o 



on 
80 


500000 


10000 


ioy 


£?no 
602 


( n 
v 


A 1 0\ 

412) 


117DZ 


no 
98 


( on 
oO 


nn\ 
00) 


A 




nn c tt c oa 
0057 (bl\) 


( n 
{{) 


AAAAAA tr 4 A 

00000054) 


4 


OA 

2y 


{ A 

v 


06) 



2 


T A 

74 


(O.Oo) 


512 


o 



yo 


500000 


10000 


1 O A 

184 


606 


/ n 
v 


444 j 


1 A A O A 


tr A 

50 


( oo 
o2 


70\ 

72) 


A 




nn c i oaTa 
0051o675 


( n 
{{) 


AAAAAA tr A\ 

00000050) 


4 


rr 
52 


( A 

v 


07) 



2 


TO 

76 


(O.Oo) 


512 


A 

4 


AA 

00 


1000000 


10000 


OAO 

208 


ono 

29o 


/ n 

(0 


o o o\ 

ooo) 


1 AAAA 

16999 


O A 

84 


/ oo 

2o 


T A \ 

79) 





A A A C TT /I 1 

00457741 


(0 


AAAAAAOOA 

OOOOOO00) 


A 

4 


66 


/ A 
( U 


05) 


2 


a a 

69 


/a ao\ 

(0.02) 


512 


A 

4 


1 a 
10 


500000 


10000 


ooo 
228 


tr A O 

54o 


/ n 
v 


a no \ 

4yo) 


1 m ot 

lyio/ 


no 

yo 


( oo 
OO 


65) 


A 




nn a ao/^ i o 
00408612 


( n 
{{) 


AAAAAA A A\ 

00000044 ) 


4 


Td 

76 


( A 

v 


07) 


2 


TA 

70 


l (\ AO\ 

(O.Oo) 


512 


4 


1 rr 

15 


200000 


10000 


ooo 

238 


TOT 

737 


/ n 

(0 


ono\ 

802) 


om i a 

20119 


A A 

40 


/" tr o 

53 


26) 





aao o rooo 

00386233 


(0 


AAAAAA/ 1 O A 

00000068) 


A 

4 


Td 

76 


/ A 

( 


12) 


2 


z?o 

62 


/a a rr \ 

(0.05) 


512 


A 

4 


20 


500000 


10000 


O A T 

247 


z^oa 

620 


/ n 

(0 


rr n 1 \ 

501) 


01 A A T 

21047 


86 


/ OO 

32 


on\ 

20) 





AAO Z" 1 tr 1 A /I 

00365104 


/A 

(0 


AAAAAA A A \ 

00000041) 


A 

4 


TA 

70 


/ A 

( 


07) 


2 


TA 

70 


/A AO\ 

(0.03) 


512 


A 

4 


OA 

30 


1000000 


10000 


0/"*0 

262 


rAA 

509 


/ n 

(0 


o rr a\ 

359) 


OO C AA 

22500 


1 A 

19 


21 


a a \ 

64) 





Aoooro tr i 

00326851 


(0 


AAAAAAOT\ 

00000027) 


A 

4 


z?o 

62 


/ A 

( 


05) 


2 


Z? O 

63 


/A AO\ 

(0.02) 


512 


A 

4 


40 


1000000 


10000 


OT A 

274 


f\ A T 

947 


/ n 

(0 


o rr a\ 

359) 


23690 


60 


/ OA 

20 


o o\ 

38) 





AAOAOO 1 O 

002y281o 


(0 


AAAAAAO tr \ 

00000025) 


A 

4 


rr 1 

51 


/ A 

( 


05) 


2 


z?o 

62 


/A AO\ 

(0.02) 


512 


A 

4 


C A 

50 


1000000 


10000 


OO A 

284 


a a n 

440 


/ n 

(0 


OZ? 0\ 

363) 


o a a 1 o 

24612 


a a 

46 


/ 1 A 

19 


rr n\ 

59) 





AAO/^O/^OO 

00262622 


/"A 

(0 


AAAAAAO A \ 

00000024) 


A 

4 


rr 

52 


/ A 

( 


05) 


2 


rr 

58 


/A AO\ 

(0.02) 


512 


A 

4 


Z? A 

60 


1000000 


10000 


OAO 

2y2 


con 

529 


/" A 

(0 


OfOA 

od2) 


O C OTT 

25377 


A A 

99 


/ 1 o 

18 


rr t\ 

57) 


u 


AAOO tT TOT 

00235787 


(0 


AAAAAAOO\ 

00000022) 


A 

4 


A A 

40 


/ A 
( U 


05) 


2 


rr z? 

56 


/A AO\ 

(0.02) 


512 


A 

4 


DO 


1000000 


10000 


2yt> 


ODl 


/ n 
v 


o^o\ 

od2) 


O C TOO 

25 f 22 


1 o 
12 


/ 1 o 

18 


on\ 

2y) 


A 




AAOOO A O 4 

0022o424 


/A 

(0 


AAAAAAOO\ 

00000022 ) 


A 

4 


/i 
42 


( A 

v 


05) 


2 


rr a 

5y 


( (\ AO\ 

(0.02) 


512 


A 

4 


7L) 


1000000 


10000 


ono 

2y& 


ac\ a 

by4 


/ n 
v 


o^n\ 

odO) 


O CAAO 

25ybo 


60 


/ IT 

17 


yi) 


A 




AAOI 1 OOO 

00211828 


Cn 
(0 


AAAAAAOI \ 

00000021 ) 


A 

4 


oa 
06 


( A 

v 


05) 


2 


rr a 
56 


/A AO\ 

(0.02) 


512 


A 

4 


on 
80 


1000000 


10000 


ono 
oOo 


oz;o 
862 


( n 
v 


o^o\ 

od2) 


ot? /i t?n 
264b0 


o c 
25 


( AT 

17 


rr j \ 
54) 


A 




AA1 A A 4 0£? 

00iy04o6 


Cn 
(0 


AAAAAAO A\ 

00000020) 


A 

4 


/l A 

40 


( A 

v 


05) 



2 


rr 
58 


( (\ AO\ 

(0.02) 


512 


A 
4t 


on 


1000000 


10000 


OUo 


00& 
ZZO 


fn 


ODO j 


ZDo i 1 


1 A 


( 17 


uy ) 


n 
u 


nm 7i ^97 

UUl 1 10Z f 


(C\ 


nnnnnni o\ 
uuuuuuiy j 


A 


OD 


y u 


05) 



z 


OD 


/n Ao^ 
(U.UZ ) 


512 


rr 



uu 


1000000 


10000 


oil 


once 

2yo 


f n 
v U 


ooy j 


071 C7 

2 / lb / 


U8 


( 1ft 

lo 


a o^ 

4z ) 


A 




nni c /i i ni 
UU154iyi 


( c\ 


0U0U0U18 ) 


4 


OA 

Z4 


v 


04) 



z 


K K 

55 


/ a no~\ 
(U.U2 ) 


512 


rr 



1 n 
10 


1000000 


10000 


olo 




f (\ 
V) 


oD2 J 


2 /40 / 


TA 
( 4 


( 1ft 

lo 


Ko^ 
52 ) 


A 




nni QQOQQ 

uulo8888 




uuouuui ( ) 


A 

4 


OA 

oy 


v 


05) 




z 


CO 

58 


/ a no~\ 
(U.U2 ) 


512 


rr 



on 
ZU 


1000000 


10000 


Q 1 ft 

olO 


1/11 
141 


1 CI 

v U 


Q^1 "\ 

oDl ) 


Z ( DoO 


10 


( 1ft 
10 


1 7\ 
1 / ) 


A 

u 


nni OKI ftQ 
UUlZOlOo 


( n 
\U 


UUUUUU10 J 


A 

4 


07 
/ 


U 


05) 




z 


KQ 

Oo 


/ a no~\ 
(U.UZ ) 


512 


rr 



on 
oU 


1000000 


10000 


ol / 


4Dy 


f n 
v U 


oDZ ) 


OTono 
Z / oUo 


A Q 
4o 


( 1ft 
10 


1 ft^ 

10 ) 


A 

u 


nni 1 oq 1 k 
UUllZolO 


f(\ 
[\J 


AAAAAA1 K ^ 

UUUUUU10 J 


A 

4 


A A 

44 


{ n 

v U 


05) 




z 


00 


/ a no~\ 
(U.UZ ) 


512 


rr 



a n 
4U 


1000000 


10000 


qi q 
olo 


K1 c 
Olo 


v U 


Q^/i ^ 

oD4 ) 


z /yoi 


ft 1 
Ol 


( 1ft 
10 


Uo ) 


A 

u 


nni ni ft07 
uuiuioy ( 




AAAAAA1 A \ 

UUUUUU14 J 


A 

4 


/l K 

40 


{ n 

v U 


05) 




z 


KQ 

Oo 


/ a no~\ 
(U.UZ ) 


512 


rr 



rcn 
5U 


1000000 


10000 


oon 
ozU 


11/1 
114 


( n 
v U 


oD4 ) 


oonoo 
zo0o2 


1 n 
10 


15 


88 ) 


A 




nnnm 701 

uouyi ( 21 




AAAAAA1 A \ 

U00UUU14 ) 


A 

4 


A A 

44 


y 


05) 




z 


CO 

5o 


/ a no~\ 
(U.U2 ) 


512 


rr 



ftn 

6U 


600000 


10000 


ozU 


/ 04 


f n 
v U 


a 7n^ 
4/0) 


OQ1 

2olDD 


Q/l 

84 


^ on 
2U 


a Q^ 
4o ) 


A 




U0082 (\)o 




UOUOUUl ( ) 


A 

4 


/i 

4z 


v 


06) 




z 


4y 


/ A C\0\ 

(U.Uo ) 


512 


rr 



/ u 


1000000 


10000 


oil 


Q A Q 

o4o 


f n 
v U 


ODO ) 


OOO /| T 

Zo24 / 


ft /i 
04 


/ iK 
15 


ft /i ^ 

64 ) 


A 




aaat/1 7/1 n 
UUU ( 4 / 4U 




UUUUUU12 ) 


A 

4 


OA 

oy 





05) 




z 


K 1 

51 


/ a no~\ 
(U.U2 ) 


1024 


o 



tr a 
50 


200000 


10000 


no 

yo 


782 


( 1 
v l 


1 ^n\ 

ioy) 


5025 


oo 
28 


/ 1 

21 


nn\ 
00) 


A 




AA O OA/" O 4 

008o06o4 


( n 
{{) 


AAAAAA tr tr A 

00000055 ) 


1 


TA 

7y 


/ A 

v 


03) 


2 


AO 

y2 


(0.06) 


1024 


A 

4 


AA 

L)L) 


200000 


10000 


O tr O 

25o 


O rr tr 
OOO 


/ 1 
v l 


iDl) 


2boob 


1 n 

iy 


115 


a a\ 
46) 


A 




A A 4 K O 1 TA 

004581 ([) 


( n 
{{) 


AAAAAAO O A 
000000O8 ) 


3 


tr 
5o 


{ A 

v 


08) 


2 


Td 

76 


( (\ a rr \ 

(0.05) 


1024 


4 


OA 

20 


200000 


10000 


351 


no o 

988 


v l 


/i O /i \ 

424) 


/i O 1 OO 

4olo2 


TC\ 

70 


l A OT 

167 


96) 





AAO/ 1 tr A A A 

00365419 


(0 


C\ A A A A A O O \ 

00000033) 


A 

4 


66 


/ A 

( 


11) 


2 


64 


/a a rr \ 

(0.05) 


1024 


4 


C A 

50 


200000 


10000 


A O 1 

481 


nn 1 

901 


/ 1 
v l 


rr oo\ 

582) 


A AAA 

64990 


c o 

58 


I A OT 

167 


on\ 

30) 





AAO/ 1 OT1 O 

00262718 


(0 


AAAAAAOT\ 

00000027) 


4 


TO 

72 


/ A 

( 


12) 


2 


rr 

58 


/a a rr \ 

(0.05) 


1024 


4 


TA 

70 


200000 


10000 


E A O 

543 


1 OO 

182 


v l 


c /i rr \ 

645) 


T A A O 1 

74421 


OT 

27 


152 


oo\ 

23) 





AAOI 1 Of 1 A 

00211861 


(0 


AAAAAAO A \ 

00000024) 


4 


66 


/ A 

( 


11) 


2 


rr t 

57 


/a a rr \ 

(0.05) 


1024 


A 

4 


on 
80 


435000 


10000 


564 


506 


( 1 
v l 


102J 


77f!7n 

77670 


7A 

79 


y5 


63j 


A 
U 


AA1 AAtT AO 

00iy0502 


(U 


AAAAAA1 C \ 
00000015 J 


4 


45 


( A 

k u 


07) 


2 


56 


/A AO\ 

(0.03J 


1024 


A 

4 


90 


350000 


10000 


581 


OO 1 

ozl 


/ 1 
v l 


233j 


onno rr 

80085 


O 1 

31 


/i no 
10o 


O A \ 

34j 


A 
U 


AA1 TA O K A 

001 7 lo5y 


(U 


AAAAAA1 Ci \ 
00000016 J 


4 


/I A 

49 


( A 

k u 


08) 


2 


57 


/A A yf \ 
(0.04J 


1024 


rr 



AA 

00 


200000 


10000 


597 


C\T A 

074 


/ 1 
v l 


646J 


oo /i m 
82401 


OA 

80 


/i oo 
128 


63j 


A 
U 


AA1 tr 4 OO/? 
00154206 


(U 


AAAAAA1A\ 
00000020J 


A 

4 


25 


( A 

k u 


10) 


2 


4/ 


/A A A \ 
(0.04J 


1024 





1 n 

1U 


500000 


10000 


OU4 


OOO 


1 




oo i uy 


AA 


( 70 


u t ) 


u 


nni ^&go^ 
uuiooyzo 


^A 


AAAAAA1 0\ 
UUUUUU1Z j 


/I 




( A 

v U 


07) 



z 


oz 


i'a n"?^ 


1024 


5 


20 


315000 


10000 


614 


282 


(i 


298) 


85148 


87 


( 97 


33) 





00125153 


(0 


00000014) 


4 


33 


(0 


08) 


2 


55 


/A A/1 \ 

(U.U4 j 


1024 


5 


30 


300000 


10000 


619 


635 


(i 


346) 


86150 


77 


( 99 


78) 





00112806 


(0 


00000014) 


4 


46 


(0 


09) 


2 


44 


/A A/] \ 

^U.U4 J 


1024 


5 


40 


500000 


10000 


624 


864 


(i 


029) 


86978 


99 


( 75 


68) 





00101735 


(0 


00000010) 


4 


45 


(0 


07) 


2 


55 


(U.Uo j 


1024 


5 


50 


500000 


10000 


630 


919 


(i 


022) 


87719 


01 


( 72 


87) 





00091759 


(0 


00000010) 


4 


32 


(0 


06) 


2 


51 


/A A 1 ^ 

(U.Uo j 


1024 


5 


60 


500000 


10000 


634 


020 


(i 


030) 


88261 


17 


( 72 


93) 





00082786 


(0 


00000009) 


4 


40 


(0 


07) 


2 


53 


(0 03^ 


1024 


5 


70 


350000 


10000 


636 


439 


(i 


243) 


88750 


82 


( 86 


74) 





00074709 


(0 


00000010) 


4 


39 


(0 


08) 


2 


53 


(0.04) 


1024 


5 


80 


500000 


10000 


639 


532 


(i 


034) 


89200 


43 


( 71 


47) 





00067443 


(0 


00000008) 


4 


39 


(0 


07) 


2 


53 


(0.03) 


1024 


5 


90 


400000 


10000 


640 


955 


(i 


148) 


89486 


27 


( 79 


44) 





00060909 


(0 


00000009) 


4 


36 


(0 


07) 


2 


54 


(0.03) 


1536 


5.40 


200000 


10000 


921.169 


(2.484) 


167548.48 


(247.24) 


0.00101730 


(0.00000011) 


4.47 


(0.11) 


2.47 


(0.04) 


1536 


5.50 


300000 


10000 


929.501 


(2.010) 


169381.15 


(196.44) 


0.00091779 


(0.00000008) 


4.45 


(0.09) 


2.56 


(0.04) 


1536 


5.60 


500000 


10000 


939.605 


(1.537) 


171375.42 


(146.36) 


0.00082793 


(0.00000006) 


4.33 


(0.06) 


2.51 


(0.03) 


1536 


5 


70 


500000 


10000 


945 


918 


(1 


553) 


172479.44 


(147.38) 


0.00074711 


(0.00000006) 


4.48 


(0.07) 


2 


51 


(0.03) 


1536 


6.00 


200000 


10000 


963.094 


(2.481) 


175739.34 


(225 


37) 


0.00055010 


(0.00000008) 


4 


39 


(0.11) 


2.46 


(0.04) 



Table 1: [Third page of 3-page table] 



41 



Total Discard 



Xstagg 



E 



Hnt.Mi 



7~int,£ 



32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 



0.50 
1.00 
1.50 
2.00 
2.50 
2.60 
2.70 
2.80 
2.90 
3.00 
3.10 
3.20 
3.30 
3.40 
3.50 
3.60 
3.70 
3.80 
3.90 
4.00 
4.10 
4.20 
4.30 
4.40 
4.50 
4.60 
4.70 
4.80 
4.90 
5.00 



1000000 
1000000 
1000000 
1000000 
1000000 
200000 
200000 
200000 
200000 
1000000 
200000 
200000 
200000 
200000 
1000000 
200000 
200000 
200000 
200000 
1000000 
200000 
200000 
200000 
200000 
1000000 
200000 
200000 
200000 
200000 
1000000 



10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 



0.408 
0.713 
0.963 
1.284 
1.494 
1.597 
1.513 
1.597 
1.707 
1.694 
1.686 
1.737 
1.727 
1.744 
1.832 
1.791 
1.855 
1.870 
1.864 
1.865 
1.837 
1.890 
1.907 
1.917 
1.907 
1.906 
1.911 
1.905 
1.996 
1.960 



0.088 
0.051 
0.039 
0.029 
0.026 
0.055 
0.060 
0.056 
0.053 
0.023 
0.053 
0.052 
0.052 
0.051 
0.022 
0.052 
0.050 
0.049 
0.049 
0.022 
0.050 
0.049 
0.048 
0.048 
0.021 
0.049 
0.049 
0.049 
0.048 
0.021 



1.64 
2.62 
3.86 
5.22 
6.45 
6.72 
6.85 
7.06 
7.32 
7.48 
7.62 
7.76 
7.90 
8.03 
8.20 
8.27 
8.44 
8.52 
8.58 
8.68 
8.69 
8.79 
8.87 
8.96 
8.99 
9.05 
9.06 
9.11 
9.20 
9.21 



0.00 
0.00 
0.01 
0.01 
0.01 
0.03 
0.03 
0.03 
0.03 
0.01 
0.03 
0.03 
0.03 
0.03 
0.02 
0.04 
0.04 
0.04 
0.04 
0.02 
0.04 
0.04 
0.04 
0.04 
0.02 
0.04 
0.04 
0.04 
0.04 
0.02 



0.33333424 
0.20768309 
0.12339135 
0.07216655 
0.04227138 
0.03809065 
0.03421029 
0.03081867 
0.02778558 
0.02498434 
0.02251044 
0.02031029 
0.01829279 
0.01646475 
0.01488912 
0.01341060 
0.01207502 
0.01094776 
0.00985439 
0.00891195 
0.00805595 
0.00726329 
0.00656180 
0.00591430 
0.00536798 
0.00486859 
0.00439399 
0.00395714 
0.00358136 
0.00323220 



0.00004069) 
0.00003587) 
0.00002967) 
0.00002316) 
0.00001777) 
0.00003844) 
0.00003588) 
0.00003480) 
0.00003256) 
0.00001345) 
0.00002917) 
0.00002774) 
0.00002590) 
0.00002463) 
0.00001037) 
0.00002191) 
0.00002158) 
0.00002028) 
0.00001904) 
0.00000792) 
0.00001717) 
0.00001631) 
0.00001557) 
0.00001478) 
0.00000617) 
0.00001338) 
0.00001268) 
0.00001187) 
0.00001127) 
0.00000472) 



2.47 
2.49 
2.54 
2.54 
2.56 
2.52 
2.72 
2.53 
2.53 
2.59 
2.56 
2.54 
2.60 
2.54 
2.63 
2.71 
2.61 
2.59 
2.63 
2.63 
2.58 
2.69 
2.57 
2.61 
2.62 
2.72 
2.62 
2.61 
2.78 
2.63 



0.02) 
0.02) 
0.02) 
0.02) 
0.02) 
0.05) 
0.05) 
0.05) 
0.05) 
0.02) 
0.05) 
0.05) 
0.05) 
0.05) 
0.02) 
0.05) 
0.05) 
0.05) 
0.05) 
0.02) 
0.05) 
0.05) 
0.05) 
0.05) 
0.02) 
0.05) 
0.05) 
0.05) 
0.05) 
0.02) 



2.86 
3.12 
3.39 
3.51 
3.55 
3.57 
3.49 
3.63 
3.51 
3.53 
3.50 
3.53 
3.45 
3.49 
3.55 
3.43 
3.64 
3.56 
3.53 
3.52 
3.53 
3.51 
3.54 
3.54 
3.54 
3.54 
3.52 
3.45 
3.43 
3.47 



64 
64 
64 
64 
64 
64 



2.50 
3.00 
3.50 
4.00 
4.50 
5.00 



1000000 
1000000 
1000000 
1000000 
1000000 
1000000 



10000 
10000 
10000 
10000 
10000 
10000 



1.559 
1.765 
1.729 
1.754 
1.890 
1.894 



0.093 
0.083 
0.085 
0.084 
0.078 
0.078 



6.47 
7.47 
8.17 
8.65 
9.00 
9.18 



0.01 
0.01 
0.01 
0.02 
0.02 
0.02 



0.04228318 
0.02497250 
0.01488177 
0.00890631 
0.00536320 
0.00323899 



0.00000895) 
0.00000678) 
0.00000517) 
0.00000397) 
0.00000308) 
0.00000236) 



2.52 
2.51 
2.51 
2.49 
2.52 
2.52 



0.02) 
0.02) 
0.02) 
0.02) 
0.02) 
0.02) 



3.60 
3.56 
3.53 
3.53 
3.56 
3.48 



Table 2: Our Monte Carlo data for the 4-state Potts antiferro magnet. "Total" is 
the total number of WSK iterations performed; "Discard" is the number of itera- 
tions discarded for equilibration. Error bars (one standard deviation) are shown in 
parentheses. 
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Lmin 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


32 


571.70, 94 


178.98, 93 


138.80, 92 


138.62, 91 


135.95, 90 


134.00, 89 


125.94, 88 




6.08, 0.0% 


1.92, 0.0% 


1.51, 0.1% 


1.52, 0.1% 


1.51, 0.1% 


1.51, 0.1% 


1.43, 0.5% 


64 


263.84, 86 


113.41, 85 


100.21, 84 


100.11, 83 


98.52, 82 


96.49, 81 


92.60, 80 




3.07, 0.0% 


1.33, 2.1% 


1.19, 11.0% 


1.21, 9.7% 


1.20, 10.3% 


1.19, 11.5% 


1.16, 15.9% 


128 


206.52, 68 


87.32, 67 


75.41, 66 


75.21, 65 


72.97, 64 


72.58, 63 


70.02, 62 




3.04, 0.0% 


1.30, 4.8% 


1.14, 20.0% 


1.16, 18.1% 


1.14, 20.7% 


1.15, 19.1% 


1.13, 22.6% 


256 


119.94, 40 


63.84, 39 


58.22, 38 


58.17, 37 


57.01, 36 


55.12, 35 


52.71, 34 




3.00, 0.0% 


1.64, 0.7% 


1.53, 1.9% 


1.57, 1.5% 


1.58, 1.4% 


1.57, 1.7% 


1.55, 2.1% 



Table 3: X 2 for the fit (|2^ ) of £(J3,2L)/£(J3,L) versus £(J3,L)/L. First line is X 2 
followed by DF (number of degrees of freedom). Second line is x 2 /F>F followed by 
the confidence level. In all cases £ min = 10. 
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p 


Lmin 




TZ for 


Xstagg,oo 


TZ for Xstagg,oo 


2.50 
2.50 
2.50 
2.50 


32 
64 

256 


13.30 (0.03) 

13.31 (0.04) 

13.32 (0.13) 
13.10 (0.42) 


0.45 (3 DF, 93.0%) 
0.27 (2 DF, 87.4%) 
0.27 (1 DF, 60.2%) 
0.00 (0 DF, 100.0%) 


200.85 (0.31) 
201.13 (0.39) 
201.51 (0.57) 
201.34 (0.80) 


2.14 (3 DF, 54.4%) 
0.84 (2 DF, 65.9%) 
0.09 (1 DF, 76.6%) 
0.00 (0 DF, 100.0%) 


2.60 
2.60 
2.60 
2.60 


32 
64 
128 
256 


15.99 (0.08) 
16.00 (0.08) 
16.00 (0.12) 
16.18 (0.36) 


0.32 (2 DF, 85.3%) 
0.31 (2 DF, 85.8%) 
0.31 (1 DF, 57.7%) 
0.00 (0 DF, 100.0%) 


272.17 (0.69) 
272.22 (0.70) 
272.54 (0.79) 
273.13 (1.10) 


1.54 (2 DF, 46.3%) 
1.38 (2 DF, 50.3%) 
0.66 (1 DF, 41.7%) 
0.00 (0 DF, 100.0%) 


2.70 
2.70 
2.70 
2.70 


32 
64 
128 
256 


19.29 (0.05) 
19.29 (0.06) 
19.27 (0.12) 
19.13 (0.33) 


0.43 (3 DF, 93.3%) 
0.27 (2 DF, 87.6%) 
0.22 (IDF, 63.8%) 
0.00 (0 DF, 100.0%) 


369.78 (0.69) 
369.51 (0.75) 
369.28 (1.09) 
369.05 (1.51) 


0.73 (3 DF, 86.6%) 
0.13 (2 DF, 93.7%) 
0.05 (1 DF, 83.0%) 
0.00 (0 DF, 100.0%) 


2.80 
2.80 
2.80 
2.80 


32 
64 
128 
256 


23.27 (0.09) 

23.28 (0.09) 

23.29 (0.12) 
23.13 (0.27) 


0.91 (3 DF, 82.2%) 
1.00 (3 DF, 80.0%) 
0.86 (2 DF, 65.2%) 
0.36 (1 DF, 55.1%) 


504.73 (0.92) 

504.74 (0.91) 
504-80 (0.98) 
504.68 (1.07) 


1.94 (3 DF, 58.5%) 
1.97 (3 DF, 57.8%) 
1.99 (2 DF, 37.0%) 

1.95 (1 DF, 16.3%) 


2.90 
2.90 
2.90 
2.90 


32 
64 

256 


28.11 (0.11) 

28.13 (0.12) 

28.14 (0.14) 

28.02 (0.24) 


1.29 (3 DF, 73.2%) 
1.34 (3 DF, 72.0%) 
0.94 (2 DF, 62.6%) 
0.59 (1 DF, 44.2%) 


690.90 (1.31) 
690.94 (1.33) 
691.27 (I.4O) 
691.04 (1.51) 


3.05 (3 DF, 38.4%) 
3.32 (3 DF, 34.4%) 
1.24 (2 DF, 53.9%) 
0.99 (1 DF, 32.0%) 


2.95 
2.95 
2.95 


32 
64 
128 


30.84 (0.18) 
30.87 (0.19) 
30.84 ( 0.20) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


804.48 (4.06) 
804.91 (4.20) 
80 4.34 {4. 35) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


3.00 
3.00 
3.00 
3.00 


32 
64 
128 
256 


34.19 (0.10) 
34.23 (0.11) 
34.16 (0.16) 
34.11 (0.23) 


2.86 (4 DF, 58.1%) 
2.89 (3 DF, 41.0%) 
2.72 (2 DF, 25.7%) 
2.51 (1 DF, 11.3%) 


952.35 (1.70) 
952.09 (1.80) 
951.17 (1.93) 
950.81 (2.12) 


3.81 (4 DF, 43.2%) 
3.45 (3 DF, 32.8%) 
2.53 (2 DF, 28.3%) 
2.08 (1 DF, 14.9%) 


3.10 
3.10 
3.10 
3.10 


32 
64 
128 
256 


41.56 (0.17) 
41.59 (0.17) 

41.57 (0.18) 
41.77 (0.25) 


1.39 (3 DF, 70.8%) 
1.23 (3 DF, 74.6%) 
1.46 (2 DF, 48.1%) 
0.01 (1 DF, 92.5%) 


1319.45 (3.44) 
1319.40 (3.43) 
1318.96 (3.57) 
1322.89 (4.06) 


3.11 (3 DF, 37.5%) 
3.36 (3 DF, 34.0%) 
4.82 (2 DF, 9.0%) 
0.48 (1 DF, 49.0%) 


3.20 
3.20 
3.20 
3.20 


32 
64 

256 


50.63 (0.16) 
50.81 (0.17) 
50.70 (0.22) 
50.72 (0.27) 


5.25 (4 DF, 26.2%) 
0.60 (3 DF, 89.7%) 
0.23 (2 DF, 89.2%) 
0.17 (1 DF, 68.3%) 


1831.78 (4.49) 
1834.33 (4.75) 
1831.27 (5.25) 
1831.56 (6.16) 


5.95 (4 DF, 20.3%) 
1.35 (3 DF, 71.7%) 
0.16 (2 DF, 92.2%) 
0.07 (1 DF, 79.5%) 


3.30 
3.30 
3.30 
3.30 


32 
64 
128 
256 


61.96 (0.22) 
62.01 (0.22) 
61.98 (0.23) 
61.89 (0.27) 


0.64 (3 DF, 88.8%) 
0.95 (3 DF, 81.4%) 
0.60 (2 DF, 74.1%) 
0.10 (1 DF, 75.8%) 


2556.28 (5.94) 
2556.71 (5.95) 
2556.23 (6.16) 
2554.73 (6.52) 


0.91 (3 DF, 82.3%) 
1.15 (3 DF, 76.5%) 
0.69 (2 DF, 70.7%) 
0.06 (1 DF, 79.9%) 



Table 4: Estimated infinite-volume correlation lengths ^ and staggered susceptibil- 
ities Xstagg,oo as a function of /?, from extrapolations using various L min . Error bars 
are one standard deviation (statistical errors only). All extrapolations use n = 5 for 
£ and n = 6 for Xstagg- indicates the residual sum-of-squares (|2.25|) for combining 
estimates from different L at the same (3; the number of degrees of freedom (DF) 
and the confidence level are indicated. Our preferred fit is shown in italics; a more 
conservative good fit is shown in sans-serif; bad fits are shown in roman. [First page 
of 4-page table] 44 



p 


Lmin 


Coo 


TZ for 


X.stagg,oc 


TZ for Xstagg,oo 


3.40 
3.40 
3.40 
3.40 


32 
64 

256 


75.63 (0.26) 
75.68 (0.26) 

15.64 (0.27) 
75.60 (0.30) 


0.05 (3 DF, 99.7%) 
0.04 (3 DF, 99.8%) 
0.01 (2 DF, 99.6%) 
0.11 (1 DF, 73.5%) 


3574.64 ( 8.59) 
3575.04 ( 8.50) 
3573.43 ( 8.99) 
3571.90 ( 9.54) 


0.29 (3 DF, 96.2%) 
0.24 (3 DF, 97.0%) 
0.02 (2 DF, 99.1%) 
0.04 (1 DF, 84.6%) 


3.45 
3.45 
3.45 


32 
64 

i25 


83.56 (0.76) 
83.75 (0.76) 
83.64 (0-78) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


4220.00 ( 47.16) 
4223.68 ( 47.21) 
4215.80 ( 47.28) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


3.50 
3.50 
3.50 
3.50 


32 
64 
128 
256 


92.94 (0.25) 
93.05 (0.26) 
93.02 (0.29) 
92.90 (0.34) 


4.88 (5 DF, 43.1%) 
1.45 (4 DF, 83.5%) 
0.68 (3 DF, 87.9%) 
0.63 (2 DF, 72.8%) 


5035.19 ( 9.43) 
5034.94 ( 9.86) 
5030.37 ( 10.68) 
5025.27 ( 12.41) 


6.50 (5 DF, 26.0%) 
1.85 (4 DF, 76.3%) 
0.37 (3 DF, 94.7%) 
0.01 (2 DF, 99.3%) 


3.60 
3.60 
3.60 
3.60 


32 
64 
128 
256 


113.85 (0.41) 
113.96 (0.42) 
113.81 (O.46) 
113.70 (0.52) 


1.82 (3 DF, 61.0%) 
1.27 (3 DF, 73.7%) 
0.08 (2 DF, 96.0%) 
0.01 (1 DF, 91.5%) 


7054.43 ( 22.62) 
7058.81 ( 23.51) 
7053.38 ( 25.92) 
7048.81 ( 30.78) 


2.11 (3 DF, 54.9%) 
1.34 (3 DF, 72.0%) 
0.02 (2 DF, 98.9%) 
0.00 (1 DF, 97.9%) 


3.65 
3.65 
3.65 
3.65 


32 
64 
128 
256 


125.80 (0.55) 
125.92 (0.57) 
125.78 (0.63) 
125.66 (0.70) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


8350.93 ( 30.21) 
8354.88 ( 31.24) 
8348.35 ( 34.51) 
8341.72 ( 40.35) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


3.70 
3.70 
3.70 
3.70 


32 
64 
i£S 
256 


139.75 (0.58) 
139.91 (0.59) 

139.76 ( 0.67) 
139.51 (0.77) 


1.07 (3 DF, 78.4%) 
1.44 (3 DF, 69.7%) 
0.64 (2 DF, 72.5%) 
0.05 (1 DF, 82.7%) 


9952.36 ( 34.98) 
9952.21 ( 35.13) 
9940.95 ( 38.57) 
9927.15 ( 45.42) 


1.32 (3 DF, 72.5%) 
1.69 (3 DF, 63.8%) 
0.82 (2 DF, 66.5%) 
0.12 (1 DF, 72.4%) 


3.80 
3.80 
3.80 
3.80 


32 
64 
128 
256 


171.68 (0.70) 
171.85 (0.73) 
171.59 (0.79) 
171.23 (0.90) 


4.96 (3 DF, 17.5%) 
3.22 (3 DF, 35.9%) 
0.18 (2DF, 91.5%) 
0.00 (1 DF, 97.6%) 


14021.81 ( 49.40) 
14009.84 ( 51.53) 
13982.10 ( 54. 81) 
13953.90 ( 61.64) 


5.11 (3 DF, 16.4%) 
3.37 (3 DF, 33.8%) 
0.27 (2 DF, 87.4%) 
0.03 (1 DF, 86.8%) 


3.90 
3.90 
3.90 
3.90 


32 
64 
128 
256 


211.00 (0.88) 
211.19 (0.91) 
210.98 (0.99) 
210.74 (1.13) 


0.98 (3 DF, 80.6%) 
1.09 (3 DF, 78.0%) 
1.30 (2 DF, 52.3%) 
0.34 (1 DF, 55.9%) 


19803.32 ( 77.69) 
19790.24 ( 81.81) 
19765.50 ( 87.33) 
19737.60 ( 99.20) 


1.08 (3 DF, 78.3%) 
1.26 (3 DF, 73.8%) 
I.48 (2 DF, 47.8%) 
0.34 (1 DF, 56.1%) 


3.95 
3.95 
3.95 


32 
64 

128 


237.95 (6.17) 
236.37 (5.97) 
235.80 (6.14) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


24159.76 (969.01) 
23848.23 (931.88) 
23773.60 ( 954.62) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


4.00 
4.00 
4.00 
4.00 


32 
64 

256 


259.76 (0.92) 
260.24 (0.96) 
259.98 (1.07) 
259.73 (1.30) 


4.79 (5 DF, 44.2%) 
1.09 (4 DF, 89.7%) 
0.92 (3 DF, 82.1%) 
0.40 (2 DF, 81.7%) 


27975.61 ( 95.23) 
27993.91 (103.35) 
27970.00 (116.37) 
27941.80 (143.30) 


5.05 (5 DF, 41.0%) 
1.27 (4 DF, 86.6%) 
1.12 (3 DF, 77.2%) 
0.58 (2 DF, 74.8%) 


4.10 
4.10 
4.10 
4.10 


32 
64 
128 
256 


317.61 (1.78) 
318.18 (1.82) 
317.70 (1.97) 
317.47 (2.39) 


0.61 (2 DF, 73.6%) 
0.28 (2 DF, 87.0%) 
0.30 (2 DF, 86.0%) 
0.16 (1 DF, 69.3%) 


39274.02 (257.27) 
39291.96 (262.01) 
39219.30 (273.74) 
39222.30 (325.32) 


0.61 (2 DF, 73.7%) 
0.30 (2 DF, 86.0%) 
0.34 [2 DF, 84.4%) 
0.18 (1 DF, 67.5%) 


4.15 
4.15 
4.15 
4.15 


32 
64 
128 
256 


355.81 (3.47) 
356.61 (3.49) 
356.11 (3.58) 
355.96 (3.93) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


47353.74 (572.90) 
47389.86 (581.04) 
47294.50 (588.26) 
47319.30 (627.61) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
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Lmin 


Coo 


71 for 


X.stagg,oo 


TZ for Xstagg,oo 


4.20 
4.20 
4-20 
4.20 


32 
64 

i2£ 

256 


394.98 ( 1.92) 
395.38 ( 1.96) 
394.85 ( g.ljg) 
394.00 ( 2.47) 


1.43 (3 DF, 69.8%) 

1.44 (3 DF, 69.7%) 
1.23 (3 DF, 7^.5%) 
1.71 (2 DF, 42.6%) 


56299.87 ( 272.76) 
56273.48 ( 284.20) 
56161.10 ( 300.19) 
56045.30 ( 332.62) 


2.48 (3 DF, 47.9%) 
2.48 (3 DF, 47.9%) 
2.26 (3 DF, 52.1%) 
3.02 (2 DF, 22.1%) 


4.30 
4.30 
4-30 
4.30 


32 
64 
128 
256 


489.45 ( 3.01) 
489.42 ( 3.14) 
48845 ( 3.56) 
488.41 ( 4.61) 


1.82 (2 DF, 40.2%) 
0.99 (2 DF, 61.1%) 
1.00 (2 DF, 60.8%) 
0.94 (1 DF, 33.2%) 


80567.72 ( 664.23) 
80425.79 ( 687.24) 
80248.30 ( 747.13) 
80379.20 ( 962.86) 


1.82 (2 DF, 40.3%) 
1.03 (2 DF, 59.7%) 
1.05 (2DF, 59.0%) 
1.00 (1 DF, 31.7%) 


4.40 
4.40 

440 

4.40 


32 
64 
128 
256 


608.66 ( 4.77) 
607.55 ( 4.90) 
606.22 ( 54 5) 
607.38 ( 7.01) 


0.50 (2 DF, 77.9%) 
0.68 (2 DF, 71.2%) 
0.75 (2 DF, 68.7%) 
0.13 (1 DF, 71.4%) 


116045.88 ( 1258.07) 
115438.34 ( 1297.88) 
115121.00 ( 1390.11) 
115580.00 ( 1724.71) 


0.50 (2 DF, 77.9%) 
0.69 (2 DF, 70.8%) 
0.76 (2BF, 684%) 
0.13 (1 DF, 71.5%) 


4.50 
4.50 
4-50 
4.50 


32 
64 
128 
256 


743.91 ( 4.99) 
741.86 ( 5.25) 
740.18 ( 5.90) 
737.83 ( 7.34) 


17.63 (5 DF, 0.3%) 
11.36 (4 DF, 2.3%) 
8.57 (3 DF, 3.6%) 
6.81 (2 DF, 3.3%) 


161644.96 ( 1441.16) 
160749.98 ( 1524.23) 
160184.00 ( 1641.99) 
159681.00 ( 1957.27) 


18.91 (5 DF, 0.2%) 
12.20 (4 DF, 1.6%) 
9.19 (3 DF, 2.7%) 
7.34 (2 DF, 2.5%) 


4.60 
4.60 
4-60 
4.60 


32 
64 
72£ 
256 


940.15 ( 9.85) 
932.11 (10.39) 
929.18 (H.49) 
929.73 (14.09) 


2.23 (2 DF, 32.8%) 

1.00 (2 DF, 60.6%) 

1.01 (2 DF, 604%) 
0.93 (1 DF, 33.6%) 


239226.36 ( 3809.46) 
235469.71 ( 3977.98) 
234461.00 ( 4288.18) 
235038.00 ( 5147.64) 


2.27 (2 DF, 32.2%) 
1.02 (2 DF, 60.2%) 
1.02 (2 DF, 60.1%) 
0.95 (1 DF, 33.1%) 


4.65 
4.65 
4-65 
4.65 


32 
64 
128 
256 


1049.65 (12.82) 
1036.57 (13.29) 
10334 9 (14.50) 
1034.53 (17.16) 


2.08 (1 DF, 14.9%) 
3.93 (1 DF, 4.8%) 
4.08(1 BF, 4.3%) 
4.14 (1 DF, 4.2%) 


287837.94 ( 5463.27) 
281335.40 ( 5552.70) 
280163.00 ( 5953.38) 
281019.00 ( 6811.95) 


2.16 (1 DF, 14.1%) 
4.10 (1 DF, 4.3%) 
4.24 (1 DF, 3.9%) 
4.32 (1 DF, 3.8%) 


4.70 
4.70 
4-70 
4.70 


32 
64 
128 
256 


1149.03 (11.81) 
1140.40 (12.25) 
1137.34 (13.39) 
1138.49 (16.23) 


0.39 (3 DF, 94.3%) 
0.16 (3 DF, 98.4%) 
0.20 (3 DF, 97.7%) 
0.20 (2 DF, 90.5%) 


334974.41 ( 5138.61) 
330406.55 ( 5250.77) 
329280.00 ( 5587.85) 
330203.00 ( 6544.22) 


0.40 (3 DF, 94.0%) 
0.17 (3 DF, 98.3%) 
0.21 (3 DF, 97.6%) 
0.21 (2 DF, 90.1%) 


4.80 
4.80 
4-80 
4.80 


32 
64 
128 
256 


1435.78 (15.18) 
1423.46 (16.09) 
1419.96 (17.36) 
1418.50 (20.45) 


3.99 (3 DF, 26.3%) 
2.25 (3 DF, 52.2%) 
2.00 (3 DF, 57.3%) 
0.45 (2 DF, 79.9%) 


485576.51 ( 7655.78) 
478214.99 ( 8078.91) 
476453.00 ( 8528.54) 
476657.00 ( 9743.68) 


3.86 (3 DF, 27.7%) 
2.09 (3 DF, 55.3%) 

1.87 (3 DF, 60.0%) 
0.45 (2 DF, 80.0%) 


4.90 
4.90 
4-90 
4.90 


32 
64 
128 
256 


1779.47 (24.92) 

1759.03 (25.37) 

1755.04 ( 26.90) 
1753.30 (30.59) 


3.54 (3 DF, 31.6%) 
3.76 (3 DF, 28.9%) 
3.83 (3 DF, 28.1%) 
3.72 (2 DF, 15.5%) 


693347.95 (15088.08) 
679244.94 (15148.94) 
677221.00 (15901.70) 
677150.00 (17675.40) 


3.66 (3 DF, 30.1%) 
3.90 (3 DF, 27.2%) 
3.97 (3 DF, 264%) 
3.86 (2 DF, 14.5%) 


5.00 
5.00 
5.00 
5.00 


32 
64 

256 


2218.69 (44.16) 
2198.02 (44.85) 
2189.50 (47.16) 
2202.41 (53.21) 


10.43 (5 DF, 6.4%) 
4.62 (4 DF, 32.8%) 
5.32 (3 DF, 15.0%) 
2.84 (2 DF, 24.2%) 


1001133.00 (32465.03) 
986315.31 (32381.73) 
980853.00 (33603.90) 
991980.00 (37473.90) 


10.42 (5 DF, 6.4%) 
4.75 (4 DF, 31.4%) 
5.53 (3 DF, 13.7%) 
2.95 (2 DF, 22.9%) 


5.10 
5.10 
5.10 
5.10 


32 
64 
128 
256 


2670.51 (51.77) 
2627.87 (51.80) 
2618. 49 (54.00) 
2619.75 (60.40) 


0.02 (3 DF, 99.9%) 
0.23 (3 DF, 97.3%) 
0.14 (3 DF, 98.7%) 
0.10 (2 DF, 94.9%) 


1371476.00 (42469.89) 
1332282.00 (42007.23) 
1325547.00 (43166.40) 
1328528.00 (47282.20) 


0.02 (3 DF, 99.9%) 
0.21 (3 DF, 97.6%) 
0.14 (3 DF, 98.7%) 
0.10 (2 DF, 94.9%) 
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Lmin 


Coo 


TZ for 


Xstagg,oo 


TZ for Xstagg,oo 


5.20 
5.20 
5.20 
5.20 


32 
64 
128 
256 


3488.77 ( 96.81) 
3431.88 ( 100.37) 
3416.96 ( 99.80) 
3419.80 ( 109.80) 


0.79 (3 DF, 85.2%) 
0.64 (3 DF, 88.8%) 
0.54 (3 DF, 91.0%) 
0.50 (2 DF, 78.0%) 


2133709.00 ( 96738.46) 
2070752.00 ( 98793.74) 
2058294.00 ( 97987.50) 
2064205.00 ( 106605.00) 


0.77 (3 DF, 85.6%) 
0.59 (3 DF, 89.9%) 
0.52 (3 DF, 91.4%,) 
0.49 (2 DF, 78.2%) 


5.30 
5.30 
5.30 
5.30 


32 
64 

i2£ 

256 


4134.51 ( 138.26) 
4061.19 ( 142.43) 
4038.28 ( 142.59) 
4057.82 ( 156.95) 


0.72 (2 DF, 69.9%) 
0.50 (2 DF, 77.9%) 
0.47 (2 DF, 79.0%) 
0.54 (2 DF, 76.4%) 


2835978.00 ( 156936.69) 
2747493.00 ( 159307.94) 
2725035.00 ( 158160.00) 
2749120.00 ( 173702.00) 


0.69 (2 DF, 70.7%) 
0.48 (2 DF, 78.6%) 
0.45 (2 DF, 79.7%,) 
0.52 (2 DF, 77.1%) 


5.40 
5.40 
5.40 
5.40 


32 
64 
128 
256 


5039.09 ( 148.02) 
4953.59 ( 148.53) 
4928.06 ( i52.6£) 
4941.46 ( 163.20) 


0.28 (3 DF, 96.4%) 
0.71 (3 DF, 87.1%) 
0.60 (3 DF, 89.7%) 
1.01 (3 DF, 79.8%) 


3955726.00 ( 190762.23) 
3831850.00 ( 188716.09) 
3804126.00 ( 192397.00) 
3823421.00 ( 204251.00) 


0.27 (3 DF, 96.6%) 
0.66 (3 DF, 88.2%) 
0.57 (3 DF, 90.4%,) 
0.92 (3 DF, 82.1%) 


5.50 
5.50 
5.50 
5.50 


32 
64 
128 
256 


6377.29 ( 217.63) 
6223.22 ( 221.23) 
6203.03 ( 224. 08) 
6197.48 ( 238.44) 


1.99 (2 DF, 37.0%) 
2.19 (2 DF, 33.4%) 
1.96 (2 DF, 37.5%,) 
2.28 (2 DF, 31.9%) 


5821262.00 ( 327977.26) 
5574691.00 ( 326750.54) 
5553237.00 ( 330812.00) 
5546998.00 ( 347228.00) 


1.99 (2 DF, 36.9%) 
2.15 (2 DF, 34.2%) 
1.92 (2 DF, 38.3%) 
2.20 (2 DF, 33.3%) 


5.60 
5.60 
5.60 
5.60 


32 
64 
128 
256 


7970.66 ( 301.47) 
7801.45 ( 309.62) 
7760.02 ( 311.34) 
7790.68 ( 326.89) 


0.04 (2 DF, 97.8%) 
0.02 (2 DF, 99.0%) 
0.01 (2 DF, 99.6%) 
0.05 (2 DF, 97.4%) 


8498373.00 ( 534485.16) 
8180645.00 ( 538772.90) 
8116724.00 ( 540 193. 00) 
8178031.00 ( 561400.00) 


0.04 (2 DF, 98.0%) 
0.01 (2 DF, 99.3%) 
0.01 (2 DF, 99.7%,) 
0.04 (2 DF, 97.8%) 


5.70 
5.70 
5.70 
5.70 


32 
64 
128 
256 


9664.89 ( 456.80) 
9531.97 ( 469.77) 
9460.68 ( 459.20) 
9548.22 ( 514.52) 


0.23 (2 DF, 89.2%) 
0.02 (2 DF, 98.9%) 
0.04 (2 DF, 97.9%) 
0.04 (2 DF, 97.8%) 


11718104.00 ( 928104.64) 
11427243.00 ( 940367.35) 
11295711.00 ( 915422.00) 
11476466.00 ( 1029147.00) 


0.22 (2 DF, 89.8%) 
0.02 (2 DF, 99.0%) 
0.04 (2 DF, 98.1%) 
0.04 (2 DF, 98.0%) 


5.80 
5.80 
5.80 
5.80 


32 
64 

i2£ 

256 


12659.85 (1409.69) 
12557.84 (1515.32) 
12416.40 (I475.7O) 
12649.50 (1593.68) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%,) 
0.00 (0 DF, 100.0%) 


18383870.00 ( 3478048.00) 
18099541.00 ( 3740026.00) 
17772288.00 ( 3613560.00) 
18354042.00 ( 3979249.00) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


5.90 
5.90 
5.90 
5.90 


32 
64 
128 
256 


14796.02 (2144.26) 
14844.78 (2373.04) 
14634.30 (2366.26) 
15023.40 (2645.59) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%,) 
0.00 (0 DF, 100.0%) 


23832311.00 ( 6025802.00) 
23910024.00 ( 6775206.00) 
23371893.00 ( 6702087.00) 
24447729.00 ( 7839348.00) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


6.00 
6.00 
6.00 
6.00 


32 
64 
128 
256 


25502.30 (7489.41) 
25937.75 (8904.71) 
25480.80 (8389.04) 
26391.60 ( 00) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 


58916945.00 (37012150.00) 
60454213.00 (45821760.00) 
58758630.00 (41692067.00) 
62389831.00 ( 00) 


0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
0.00 (0 DF, 100.0%) 
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Lmin 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


32 


1738.06, 94 
18.49, 0.0% 


288.25, 93 
3.10, 0.0% 


277.13, 92 
3.01, 0.0% 


244.18, 91 
2.68, 0.0% 


243.54, 90 
2.71, 0.0% 


236.75, 89 
2.66, 0.0% 


228.29, 88 
2.59, 0.0% 


64 


932.48, 86 
10.84, 0.0% 


196.30, 85 
2.31, 0.0% 


193.24, 84 
2.30, 0.0% 


178.47, 83 
2.15, 0.0% 


177.49, 82 
2.16, 0.0% 


171.16, 81 
2.11, 0.0% 


164.95, 80 
2.06, 0.0% 


128 


692.28, 68 
10.18, 0.0% 


142.15, 67 
2.12, 0.0% 


138.73, 66 
2.10, 0.0% 


131.45, 65 
2.02, 0.0% 


128.69, 64 
2.01, 0.0% 


127.09, 63 
2.02, 0.0% 


123.49, 62 
1.99, 0.0% 


256 


387.08, 40 
9.68, 0.0% 


106.87, 39 
2.74, 0.0% 


106.47, 38 
2.80, 0.0% 


101.90, 37 
2.75, 0.0% 


100.23, 36 
2.78, 0.0% 


95.69, 35 
2.73, 0.0% 


90.82, 34 
2.67, 0.0% 



Table 5: X 2 for the fit ( ggg) of Xsta 9g (/3, 2L) /XstaggW, L) versus £(J3, L)/L. First line 
is x 2 followed by DF (number of degrees of freedom). Second line is x 2 /E ) F followed 
by the confidence level. In all cases £ min = 10. 
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FSS Plot for Correlation Length 




0.0 0.2 0.4 0.6 0.8 

£(L)/L 



Figure 1: £(/3, 2L)/£(/3, L) versus £{(3,L)/L. Symbols indicate L = 32 (+), 64 
(x), 128 (□), 256 (O), 512 (O)- Error bars are one standard deviation. Curve is a 
fifth-order fit in ( |2.22|) , with £ min = 10 and L min = 128. 
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Deviations from FSS Fit for Correlation Length 
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Deviations from FSS Fit for Correlation Length 
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Figure 2: Deviation of points from fit to F% with s = 2, n = 5, £ m j n = 10 and 
L min = 128. Symbols indicate L = 32 (+), 64 (x), 128 (□), 256 (O), 512 (O)- Error 
bars are one standard deviation. Curves near zero indicate statistical error bars (± 
one standard deviation) on the function F^(x). Plot (a) shows all points; plot (b) is 
a blow-up showing the L = 32 points at x > 0.50. 
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0.00 



Deviations from FSS Fit for Susceptibility 
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Figure 3: Deviation of points from fit to F Xstagg with s = 2, n = 6, £ m j„ = 10 and 
L min = 128. Symbols indicate L = 32 (+), 64 (x), 128 (□), 256 (O), 512 (O)- Error 
bars are one standard deviation. Curves near zero indicate statistical error bars (± 
one standard deviation) on the function F Xstagg (x). Plot (a) shows all points; plot (b) 
is a blow-up showing the L = 32 points at x > 0.50; plot (c) is a blow-up showing all 
points with x > 0.58. 
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FSS Plot for Susceptibility 




Figure 4: Xsta gg {P, 2L)/x s t agg (P, L) versus £(/3, L)/L. Symbols indicate L = 32 (+), 
64 (x), 128 (□), 256 (O), 512 (O)- Error bars are one standard deviation. Curve is 
a sixth-order fit in ( |2.22j ), with £ min = 10 and L min = 128. 
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d=2, q = 3 AF Potts Correlation Length 
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Figure 5: Infinite-volume correlation length ^ divided by e 2/3 (3 p for (a) p = 0, 
p — 1/2, (c) p — 1. Error bars are one standard deviation. 
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Figure 6: £oo/ (e 2f3 j3 p ) with p — 1, plotted versus e Note the nearly linear behavior, 
in good agreement with (|T5]). Straight line is £oc/{e w (3) = 0.01814 + 0.20051e -/J , 
which is the least-squares fit to the data with j3 > 2.95 (e _/3 ^ 0.052). 
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d = 2, q=3 AF Potts Correlation Length 
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Figure 7: £oo/e 2/3 plotted versus e~ A/3 . (a) A = 1. Straight line is £oo/e 2/3 = 0.1098 - 
1.6574e _/? , which is the least-squares fit to the data with (3 > 4.50 (e^ < 0.011). (b) 
A = 0.5. Straight line is We 2/3 = 0.1225 - 0.2929e" a5/3 , which is the least-squares fit 
to the data with (3 > 4.50 (e~°- 5/3 < 0.105). 



55 



d=2, q = 3 AF Potts Correlation Length 
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d = 2, q = 3 AF Potts Correlation Length 
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Figure 8: W(e 2/3 /3 p ) plotted versus 1/(3 for (a) p = 0, (b) p = 1/2, (c) p = 1. 
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d = 2, q=3 AF Potts Correlation Length 
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(1 = 2, q = 3 AF Potts Staggered Susceptibility 
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Figure 10: Xstagg,oo/^ 3 plotted versus (3. Error bars are those given by the triangle 
inequality, reduced by a factor of 10 for visual clarity. 
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d=2, q = 3 AF Potts Staggered Susceptibility 
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d=2, q=3 AF Potts Staggered Susceptibility 
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Fi gure 11; Infinite- volume staggered susceptibility Xstagg, oo 

divided by e (w /^l3 q for 
(a) q — 0, (b) q — 5/6, (c) q = 5/3. Error bars are one standard deviation. 
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d = 2, q = 3 AF Potts Staggered Susceptibility 
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Figure 12: Xstagg, oo /{eWWp) with g = 5/3, plotted versus e The behavior 
is reasonably linear for e _/3 < 0.03, in good agreement with (|4.8|) . Straight line is 
Xsta 9S ,oo/(e (10/3)/3 /5 5/3 ) = 0.00329 + 0.06661e" /3 , which is the least-squares fit to the 
data with (3 > 3.60. 
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d = 2, q = 3 AF Potts Energy 
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Figure 13: Energy per site E divided by e 13 , plotted versus e Symbols indicate 
L = 32 (+), 64 (x), 128 (□), 256 (O), 512 (O), 1024 (*), 1536(tt). Error bars are 
invisibly small. The uppermost points at each j3 represent the infinite- volume limit. 
Note the nearly linear behavior for small e -13 , in good agreement with ( [4.9| ). Straight 
line is E = 0.21777e _/J + 1.65303e~ 2/3 , which is the least -squares fit to the data with 
L = 1024 and (3 > 5.00. 
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Figure 14: Dynamic finite-size-scaling plot of r int ^ M 2 versus £(L)/L, assuming 
dynamic critical exponent z int M 2 = 0. Symbols indicate L = 32 (+), 64 (x), 128 
(□), 256 (O), 512 (O), 1024 (*), 1536(s). Error bars are one standard deviation. 
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Dynamic FSS for T intE : Exponent z = 




0.0 0.2 0.4 0.6 0.8 

£(L)/L 



Figure 15: Dynamic finite-size-scaling plot of Ti nt £ versus £(L)/L, assuming dynamic 
critical exponent z intt s = 0. Symbols indicate L = 32 (+), 64 (x), 128 (□), 256 (O), 
512 (O), 1024 (*), 1536(hEh). Error bars are one standard deviation. 
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Figure 16: Plot of M 2 (t) versus t/r, ;rit M 2 , using all data points. Symbols 

J stagg J stagg v y ' ,y stagg 

indicate L = 32 (+), 64 (x), 128 (□), 256 (O), 512 (O), 1024 (*), 1536(ft). Error 
bars are omitted. 
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Figure 17: Plot of p M 2 M 2 (t) versus t/T irif M 2 , subdivided by ranges of 

J stagg J stagg v f if* l stagg 

i{L)/L: (a) 0.0-0.1, (b) 0.1-0.2, (c) 0.2-0.5, (d) 0.50-0.54, (e) 0.54-0.58, (f) 0.58- 
0.63. Symbols indicate L = 32 (+), 64 (x), 128 (□), 256 (O), 512 (O), 1024 (*), 
1536(hEh). Error bars are omitted. The straight lines correspond to a pure exponential 
decay r intMltagg = r expM 2 tagg . 
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d=2, q=3 AF Potts, Relative Variance-Tinie Product 
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d=2, q=3 AF Potts, Relative Variance-Time Product 
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Figure 18: Relative variance-time product [including errors of types (i) + (ii) + (iii)] 
divided by £oo(/3) 2 , plotted versus £oo(/3)A, for two-dimensional three-state Potts 
antiferromagnet. (a) is for O — £, (b) is for (9 = Xstagg- Symbols indicate L = 128 
(□), 256 (O), 512 (O), 1024 (*), 1536(*). 
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d = 2, q = 3 AF Potts model, £ min = 20 
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Figure 19: Traditional finite-size-scaling plot of £(L)/e 2/3 versus £(L)/L, using all 
data points with £(L) > 20. Symbols indicate L = 32 (+), 64 (x), 128 (□), 256 (O), 
512 (O), 1024 (*), 1536(ft). Error bars (one standard deviation) are in almost all 
cases smaller than the symbol size. 



67 



